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NUMERICAL  INTEGRATION  OF  SYSTEMS  OF  DIFFERENTIAL  EQUATIONS  OF 
CHEMICAL  KINETICS. 


I, :  Methods  of  Integration  of  equations  of  chemical  kinetics 
using  a  computer. 


As  is  well  known,  the  mathemematieal  description  of  chemi¬ 
cal  reactions  uses  the  derivatives  of  concentration  (the  rate 
of  the  reaction)  and  their  relationship  to  the  numerical  val¬ 
ues  for  the  concentrations  themselves.  The  choice  of  differ¬ 
ential  and  not  integral  equations  for  the  description  of 
chemical  processes  is  based  on  the  consideration  that  the 
quantities  Cj(t)  have  a  completely  clear  physical  meaning, 
while  the  parameters  of  the  integral  equation  establish  only 
the  fact  that  the  concentrations  themselves  are  some  sort  of 
functions  of  time:  cj  =  cj(t),  where  j  =  1,  2,...,n,  n  being 
the  number  of  components  taking  part  in  the  reaction. 

Coefficients  called  constants  of  speeds  of  chemical  reac¬ 
tions  are  a  part  of  sets  of  differential  equations  which  fix 
the  dependence  in  time  of  c*  on  other  concentrations.  With¬ 
out  the  use  of  a  computer  tne  solution  of  such  sets  of  dif¬ 
ferential  equations  for  complex  chemical  reactions  is,  as  a 
rule,  impossible.  The  differential  equations  are  replaced 
with  sets  of  algebraic  equations  using  the  method,  employed 
in  chemical  kinetics,  of  stationary  concentrations.  In  this 
method  all  the  c<’s  are  assumed  equal  to  zero.  In  doing  this 
the  unfolding  ofJthe  process  in  time  is  lost,  as  well  as  the 
possibility  for  investigating  the  influence  of  various  con¬ 
ditions  (in  particular,  the  concentration  of  various  compo¬ 
nents  in  the  initial  moment,  temperatures,  and  so  forth). 

The  development  of  effective  methods  of  numerical  inte¬ 
gration  of  sets  of  differential  equations  describing  complex 
chemical  reactions,  proceeding  with  the  final  rates,  can  per¬ 
mit  the  solution  of  these  problems.  The  discovery  of  a  solu¬ 
tion  for  the  set  of  differential  equations  of  the  kinetics 
of  complex  chemical  reactions  gives  the  possibility  of  ob- 
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taining  the  following  information; 

1.  When  the  values  of  the  constants  of  reaction  rates 
for  various  stages  of  a  complex  reaction  are  known,  it  is 
possible  in  a  reasonable  amount  of  computer  time  to  examine 
the  various  alternative  mechanisms  hypothesized  for  the  reac¬ 
tion. 

2.  For  each  of  the  mechanisms  examined,  it  is  possible 
to  obtain  the  values  for  the  concentrations  of  all  components 
as  a  function  of  time  for  the  entire  duration  of  the  reac¬ 
tion  (for  a  given  sum  total  of  constants  of  reaction  rates). 
Therefore,  the  comparison  of  calculated  and  experimentally- 
derived  values  for  concentrations  is  possible  for  any  one  of 
them  for  any  interval  of  time,  provided  that  their  measure¬ 
ment  is  experimentally  possible. 

3.  In  the  mathematical  model  corresponding  to  each 
known  mechanism  of  a  cqmplex  chemical  reaction,  it  is  possible 
to  learn  the  degree  of  sensitivity  of  that  mechanism  to  chan¬ 
ges  in  the  constants  of  reaction  rates  in  separate  stages  of 
the  reaction  and  in  their  different  groups,  and  also  to  learn 
the  mechanism's  degree  of  sensitivity  to  varied  conditions  af¬ 
fecting  the  way  the  reaction  proceeds  (temperature,  pressure, 
Initial  distribution  of  concentrations,  etc.).  Actually, 
such  research  may  be  conducted  without  any  supplementary  ex¬ 
perimental  data  at  all.  Moreover,  using  special  mathemati¬ 
cal  methods  (the  gradient  method  and  others),  if  the  experi¬ 
mental  data  is  reworked,  it  is  possible  to  find  approximate 
values  for  the  constants  of  rates  of  reaction  as  well  as  the 
minimum  number  of  experiments  needed  to  obtain  the  constants 
to  a  given  degree  of  accuracy. 

However,  as  will  be  shown  in  detail  below,  when  using  a 
computer  to  integrate  sets  of  differential  equations  of  the 
kinetics  of  complex  homogeneous  isothermal  reactions,  con¬ 
siderable  mathematical  difficulties  arise  in  the  case  —  im¬ 
portant  in  chemistry  —  when  chemically  active  particles  such 
as  atoms  and  radicals  take  part  in  the  reaction.  This  is  so 
because  the  constants  of  rates  of  reaction  for  such  particles 
are  several  orders  of  magnitude  larger  than  the  other  con¬ 
stants  of  rates  or  reaction  going  into  the  equations.  In 
this  case  the  use  of  the  usual  methods  of  numerical  inte¬ 
gration  on  a  computer  turn  out  to  be  extremely  difficult  to 
put  into  practice. 

It  is  possible  to  formulate  the  equations  of  chemical 
kinetics  for  an  isothermal  homogeneous  reaction  (in  closed 
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systems)  in  the  following  form*: 

ft  ft 

cj  ~  2  2  ajnpcmep  +  2  Pjpfp  0) 

«—!  P<m  P—1 

(<>*»,  /  =  1,2 . n). 

where  ci  are  the  dimensionless  concentrations  connected  with 
the  dimensional  multipliers  Co  =  Tq/VqT,  where  Tg  =273°  K, 

Vo  =  22.4  liters/mole,  and  T  is  the  temperature  L°Cj ;  the  co¬ 
efficients  <*jinp  and  £  jp  which  are  not  equal  to  zero,  are 
the  corresponding  constants  of  rates  of  second  and  first  or¬ 
der  reactions,  brought  to  the  dimensionality  sec.-l;  and  t0  is 
the  time  of  the  start  of  the  reaction. 


*  For  simplicity,  the  rare  zero-order  and  third-order  reac¬ 
tions  are  here  omitted.  Their  inclusion  in  equation  set  (1) 
does  not  introduce  material  changes  in  the  discussion  to  fol¬ 
low.  It  must  be  noted  that  the  form  given  above  is  the  uni¬ 
versal  and  sole  form  of  notation  of  equations  of  chemical 
kinetics.''  We  also  note  that  generalization  of  the  method 
offered  below  to  the  case  of  open  systems  (taking  into  ac¬ 
count  hydrodynamic  phenomena  and  non-isothermality)  does  not 
entail  material  difficulties. 


The  quantities  going  into  the  set  of  equations  (1)  satisfy 
the  following  conditions:  1)  Oicj  ^  l(j  =  l,2,,..,n); 

2)<*jmp  and^jp  are  symmetric  with  respect  to  permutations  of 
the  indices  m  and  p;  3)  if  even  one  of  the  indices  m.  and  £  is 
equal  to  then  tEe  corresponding  coefficients,  if  they  are 
different  from  zero,  are  negative.  Algebraic  equations  of 
balance  enable  us  to  lower  the  order  of  the  set  of  equations 
(1). 


Grouping  together  in  each  equation  of  set  (1)  the  terms 
containing  Cj,  we  shall  rewrite  the  set  in  the  following  form: 

ci(0  =  — MOCj  +aj(<)  (/  =  1.2 . ft),  (2) 

uhete. 

n  n 

•jW  =  2  2  aimpcmcp  4"  2  P;pcp  ^ />  P  /)» 

1  p<m  p— l 

n 

(*)  =  2  I  **iJP  I  eP  +  !  Pii  I- 
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It  is  evident  that  b.r(t).  as  opposed  to  a±(t),  may  depend 
linearly  on  ci(t).  Since  aut),  bi(t)  are  the  sums  of  the 
non-negative  terms,  then  ajtt^A.^O,  and  b, (t)  ^  B,  ^  0, 
where  A.  and  B4  are  constants  which  are  the  lower  boundaries 
of  the  functions  aj(t),  bj(t)  when  t  >t0. 

For  example,  for  reactions  of  the  first  order, 

bj(t) 

1.  The  content  of  the  problem. 

The  use  of  computers  to  solve  problems  of  chemical 
kinetics  has  grown  sharply  recently  (see,  for  instance, 
sources  [1-53 ) .  However,  as  has  already  been  mentioned,  the 
use  of  standard  methods  of  numerical  integration  for  sets  of 
equations  describing  complex  chemical  reactions  meets  with 
various  difficulties.  This  is  connected  with  the  fact  that 
source  1  when  reagents  with  substantially  different  reaction 
capability  take  part  in  the  reaction  (constants  of  rates  of 
reaction  which  when  converted  to  the  same  dimensionality,  dif¬ 
fer  by  ten  orders  of  magnitude  or  more) ,  certain  of  the  equa¬ 
tions  (e.g.  those  describing  the  kinetics  of  free  radicals) 
may  in  the  course  of  numerical  integration  become  unstable 
with  the  passage  of  time,* 


*This  question  will  be  examined  in  detail  in  section  1-3 
below. 


This  peculiarity  is  specific  to  these  equations  of  chemical 
kinetics  inasmuch  as,  from  the  point  of  view  of  the  theory  of 
sets  of  ordinary  differential  equations  when  Bj  *  0,  they  are 
absolutely  stable. 

Actually,  we  shall  examine  the  conduct  of  a  line  of  inte¬ 
gration  of  one  of  the  equations  of  the  form  (2) 

J-(t)  =  a(t)  -  b(t,y)v.  We  shall  define  (t)  «  y*(t)  -  y?(t), 
wnere  y. (t)  (i  -  1,2)  is  the  basic  component  of  the  solution 
of  this  equation,  corresponding  to  the  various  given  initial 
conditions.  Adding  and  calculating  b(t,y2)y2,  we  get: 

<f>  (I)  =  -  b  (t,  vi)  v-vtlb  (l,  m)  -  b  (t,  Val¬ 
inas  much  as,  in  agreement  with  the  above,  d)>(t,y)/dy  ^  0, 
then  when  yi  >  y2>  that  is  <fe>0,  we  have  $<0  and  vice 
versa.  Thus,  the  integration  series  converges  when  inte¬ 
gration  is  carried  out  from  left  to  right,  and  as  the  mag¬ 
nitude  of  b(t,y)  increases,  stability  around  the  initial  con¬ 
ditions  increases  (source  23. 


_  k  - 
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However,  as  will  be  shown  below,  when  large  enough  values 
of  b(t,y)  are  reached,  oscillations  of  the  solution  of  the 
equation  under  examination  occur  during  the  process  of  numer¬ 
ical  integration _by  usual  methods.  The  amplitude  of  the  osci¬ 
llations  depends  on  the  path  of  integration.  As  a  result  of 
this  occurrence.the  integration  of  a  set  of  such  differential 
equations  of  chemical  kinetics  by  means  of  standard  methods 
of  numerical  analysis  is  in  practice  impossible. 


Until  now  works  have  been  published  on  various  problems 
in  chemical  kinetics  in  which  the  difficulties  of  numerical 
integration  of  differential  equation  set  (1)  in  the  case  of 
complex  chemical  reactions  are  noted  sources  l£-43 .  Hence, 
source  [2J  conducted  a  qualitative  investigation  of  the  car¬ 
rying  out  of  integral  curves  for  an  equation  of  the  form: 


{/  = 


»-C(h 


(4) 


where  Q(t,y)  and  G(t)  are  the  given  functions,  and  Q(t,y)^^, 
with |  quite  small  (for  the  purposes  of  chemical  kinetics 
IQK  ^  ,  where  'X  is  the  time  of  the  reaction). 


It  is  evident  that  this  equation  is  a  particular  case 
of  a  class  of  problems  in  chemical  kinetics,  since  it  is  es¬ 
sential  here  that  Q  =  Q(t)  <  0  and  moreover,  it  excludes, 
for  Instance,  reactions  of  recombination.  In  source  2  it 
is  noted  that  when  Q  <.  0  the  integration  series  during  inte¬ 
gration  from  left  to  right  approaches  y  =  G(t),  and  the 
speed  of  approach  sharply  increases  with  the  growth  of^Ql”1, 
As  a  result  of  the  investigation  in  source  2  a  different 
method  of  solution  of  such  equations  (called  "stiff"  equa¬ 
tions  by  the  authors  of  source  £2]  )  is  suggested.  This 
method  is  interesting  for  the  theory  of  electric  circuits  and 
for  rocket  guidance  problems,  when  Q>>  0  and  the  integral 
curves,  when  y  p  G(t),  increase  sharply  in  the  vicinity  of 
the  modulus.  But  as  applied  to  the  problems  of  chemical  ki¬ 
netics,  this  method  has  not  been  developed. 


The  difficulties  which  we  have  described  in  numerical  in¬ 
tegration  of  problems  of  the  kinetics  of  complex  chemical  re¬ 
actions  led  us  to  attempts  at  approximate  solutions  on  a  com¬ 
puter  by  means  of  linearizing  the  equation  set  (1)  (source  £)]); 
such  linearization  is  usually  undertaken  with  the  aim  of  ana¬ 
lytical  investigation  of  the  solution.  Obviously,  without 
obtaining  an  estimate  of  the  accuracy  of  the  solution  when  the 
set  is  linearized,  this  method  can  not  be  used  for  numerical 
integration  by  computer. 

Attention  is  paid  in  source  to  the  peculiarities  of 
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behavior  of  set  (2);  in  source  (V]  a  method  of  computation** 
is  suggested,  without  any  mathematical  substantiation,  which 
is  in  many  ways  analogous  to  the  method  developed  in  article 
ClJ  and  similar  to  the  one  outlined  in  the  following  section 
of  this  paper. 


•Source"  UK]  became  known  to  us  after  the  method  presented  by 
us  had  already  been  worked  out  and  confirmed  on  several  exam¬ 
ples. 

**In  this  method  an  idea  is  used  which  had  already  been  for¬ 
mulated  by  Bodenstein  [53  , known  as  the  method  of  quasi^-sta- 
tionary  (or  partially  stationary)  concentrations  (source  [63). 

Bodenstein's  idea,  more  simplified,  was  used  in  Blackmore’s 
work  D.H.  Blackmore.  Trans.  Faraday  Soo..  62,  518,  415  (1966)  . 
In  this  work  the  differential  equations  of  chemical  kinetics 
for  components  with  a  high  reaction  capacity  were  immediately 
replaced  by  algebraic  equations.  Such  a  replacement  is  not 
always  correct.  It  was  shown  by  us  under  what  conditions  such 
a  replacement  may  be  done  correctly. _ 

Thus,  in  source  [4-3  ,  it  is  asserted  a  priori  that  in  the  equa¬ 
tions  for  certain  concentrations  which  are  assumed  to  be  small 
in  the  course  of  the  whole  reaction,  the  derivatives,  starting 
from  a  certain  moment  in  time,  will  also  be  quite  small.  Con¬ 
sequently,  during  the  process  of  calculation  on  the  computer 
program,  this  moment  in  time  is  determined  and  the  right-hand 
parts  of  the  corresponding  equations  are  assumed  to  be  approxi¬ 
mately  zero. 

Thus,  in  source  [4]  ,  the  reasons  for  the  peculiarities 
of  behavior  of  the  numerical  solution  of  equations  of  kinetics 
for  complex  chemical  reactions  remained  undiscovered,  which 
led  to  the  fact  that  the  authors  of  that  source  erroneously 
opened  up  the  prospect  of  solving  beforehand  (without  the 
results  of  a  numerical  integration)  those  equations  which  can 
make  the  set  unstable  during  the  process  of  integration. 

Therefore  in  source  [4]  it  was  impossible  to  determine  the 
region  of  useability  of  the  calculation  method,  or  to  give 
an  evaluation  of  its  accuracy.  Because  of  all  this,  the  use 
of  the  method  is  made  more  difficult  and  in  certain  situations, 
(for  example,  for  second-order  reactions,  see  below)  can  lead 
to  significant  mistakes. 

2.  Details  of  the  solutions  of  equation  set  (1)  when 
there  are  particles  of  different  reaction  capability  taking 
part  in  the  complex  chemical  reaction. 
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We  shall  explain  the  nature  of  the  influence  of  chemically 
active  components  on  the  solution  of  equation  set  (1);  for 
simplicity's  sake,  we  *hall  first  examine  the  case  where  one 
component — the  k-th — has  a  high  reaction  capacity  (defini¬ 
tion  of  sufficient  conditions  for  this  to  apply  are  given 
after  inequalities  (4)  and  (5)  below). 

Let  us  assume  that  all  concentrations  c*(t)(j  =  l,2,,„.n) 
for  all  times  t  >  0  were  somehow  previously ’■'determined  for  us. 
Then  the  k-th  equation  in  set  (1),  in  agreement  with  equation 
set  (2),  may  be  written  in  the  form 

Ck  (0  ~  —  bk  ( t )  ek  -f  ak  (t)'  (3) 

Here  b.  (t)  and  av(t),  determined  by  the  k-th  equation  of 
set  (1;  and  by  the  known  values  of  the  concentrations,  are 
assumed  to  be  known  functions  of  time.  Now  with  the  aid  of 
equation  (3),  which  we  set  up  only  for  the  purpose  of  quali¬ 
tative  investigation,  we  shall  obtain  an  evaluation  for  tlae 
analysis  of  ck(t). 

It  isr  evident  that  for  a  particular  choice  of  ck(t)  for 
the  equation  being  examined,  two  cases  are  possible* 

Case  1)  The  following  conditions  hold  for  the  entire  du¬ 
ration  of  the  reaction,  whose  total  time  is  yQ: 

bk  (t)  >  h  (0  </  -  1,2,...,  k  -  1,  k  +  1,...,  n),’  (4) 

Bk>^~,  (5 V 

(the  latter  of  the  inequalities  in  (5)  derives  from  the  for¬ 
mer;  £  is  the  desired  degree  of  accuracy*) 


^Because  the  values  of  the  various  concentrations  change 
widely  (over  a  range  from  0  to  unity) ,  the  magnitude  of 
must  be  chosen  quite  small  compared  with  unity. _ 

From  an  examination  of  set  (2)  together  with  conditions  (4) 
and  (5)  it  follows  that  components  participating  in  a  first- 
order  reaction,  and  having  relatively  high  constants  of  reac¬ 
tion  rate  i<£kk!  ,  are  intentionally  included  in  this  case. 

Case  2)  It  is  not  possible  to  confirm  ahead  of  time  that 
conditions  (4)  and  (5)  are  fulfilled  over  the  whole  duration 
of  the  reaction,  but  they  at  least  hold  for  certain  parts  of 
the  reaction.  In  this  casehigh  reaction  capabilities  of  com¬ 
ponents  may  occur,  for  instance,  in  reactions  of  the  second 


order.  Then  whether  or  not  conditions  (4)  and  (5)  hold  de¬ 
pends  on  the  magnitude  of  the  value  of  the  concentrations 
having  large  constants  of  rates  of  reaction- as  co-factors 
in  the  expression  for  function  b^ft). 

In  the  simplest  case,  (which  has  been  examined  in  source 
[6 J ,  for  example),  where  ( t )  =  ak(t)  =  0,  the 

solution  of  the  equation  obtained  from  equation* (3)  has  the 
form:  - 

«»(<)  =  c*(0)  oxp  (-;//»  0  +  ^  11  -  o*P  <-  Okt)  ]. 


Since  0  ^  c  (0)  4  1,  then  when  conditions  (4)  and  (5) 
are  fulfilled,  the  first  item,  which  expresses  the  influence 
of  the  initial  conditions,  may  be  excluded  from  considera¬ 
tion  when  t 
Hence  when  t  > 


p  =  that  is,  when  expf-B^t)^^ 

£  we  have  within  the  limits  of  accuracy: 


Ct(t)tS!Ct  I 


At 

Bk 


#=®. 


where  c,(t)  and  consequently,  c.  (t)  will  be  small  compared  to 
unity.  *The  concentration  ck,  Which  we  have  obtained  in  this 
particular  case  we  have  been  examining,  is  known  as  the  quasi¬ 
stationary  concentration,  and  is  widely  used  in  chemical  kine¬ 
tics. 


An  analysis  of  the  behavior  of  Cu(t)  in  the  general 
case,  where  a,  =  a.  (t)  and  b^  =  b.  (tj  are  determined  by  sol¬ 
ving  the  whole  setKof  equations  or  chemical  kinetics,  presents 
considerable  difficulties.  Hov/ever,  for  components  whose 
equations  satisfy  conditions  (4)  and  (5),  and  in  the  case  of 
equation  set  (1),  quasi -stationary  concentrations  also  occur. 
It  is  evident  that  the  latter  concentrations  will  be  substan¬ 
tially  different  from  those  examined  earlier  in  the  particu¬ 
lar  case.  In  fact,  from  equation  (3),  inasmuch  as  a^ft)  and 
bj^t)  are  assumed  to  be  known  functions  of  time, 

«*  (0  =  e*  (0)  «p  [—  $  bk  (s)  ds]  +  $  a„  (s)  exp  [—  $  bk  (r)  dr\  ds.  (6) 

•  i  i 

Hence,  similar  to  the  particular  case  examined  above, 
the  first  item,  which  directly  reflects  the  influence  of  the 
initial  conditions,  becomes  small  enough  to  ignore  when 
t  ^  t0  -  |ln£\  /B.  .  From  the  boundedness  of  a^(t),  bj.(t)  it 
follows  that  it  is  always  possible  to  also  choose  d^  =  const, 
and  g|  =  const,  (i  =  1,2)  so  that  d^  a^(t)<  d2,  and 
£i  ^  b^(t) <.  g2,  where  d^1>  0,  g.  >  0  ana,  in  accordance 
with  inequality  (5),  gj  >  Bk. 
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Then: 


'  * 

<*!  [1  —  exp  (—  git) ]  <  J  a*  (s)  exp  [—  £&*(/•) </r] ds < 


where : 


From  equation  (6),  for  t  >  t0,  in  accordance  with  the 
foregoing,  and  within  the  accuracy  limits,  we  have 

Ox  <**  (0  <<*»• 

Inasmuch  as  <%<<^l(i  =  1,2)  in  accordance  with  inequality 
(5),  we  obtain  the  result  that  the  concentrations  of  the 
highly  reactive  elements  that  we  have  been  examining  become 
small  in  relation  to  unity  and  consequently,  the  change  in 
ak(t)  and  b,(t)  for  t>  t  is  basically  determined  by  the  con¬ 
centration  of  the  substances  having  a  relatively  small  capa¬ 
city  of  reaction  compared  to  c  (t).  Then,  using  the  continu¬ 
ous  nature  of  the  solution  of  equation  (3).  a  6>0  is  delibera 
tely  chosen  such  that,  when  t  6  (t,  t  +6  )  we  have  occurring: 

|«j --£*{<)  I  <«(<  =  1,2), 


where 


**(0  = 


*Jr<0 

m- 


*If  b^  =  b^Ct.Ck),  then  <5^  is  found  in  a  similar  fashion. 


Thus,  if  in  equation  set  (l)  there  is  one  equation  satis¬ 
fying  conditions  (h)  and  (5)  over  the  duration  of  the  entire 
reaction,  then  for  all  t  |ln£.|  /g,,  the  error  made  in  substi¬ 
tuting  the  corresponding  concentration  c.  (t)  for  ^(t)  becomes 
less  than  the  given  maximum  allowable  error;  and  in  accor¬ 
dance  with  equation  (3)  St,(t)  becomes  as  close  to  zero  as 
desired.  Characteristically,  in  this  event  the  individual 
terms  in  the  right-hand  part  of  equation  (3)  are  large  com¬ 
pared  with  the  magnitude  of  CiK(t). 

From  now  on  we  shall  call  the  function  c.  (t)  the  quasi¬ 
stationary  solution  of  equation  (3).  It  Is  evident  that  the 
quasi -stationary  situation  for  individual  substances  that  are 
participating  in  the  reaction  is  typical  for  reactions  where 
there  are  stages  with  constants  of  reaction  rates  sufficiently 
large  and  varied  in  magnitude.  Al'co,  the  quasi-s  tat  ionary 
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values  change  over  time  as  a  function  of  the  behavior  of  the 
separate,  less  reactive  components  of  reaction. 

In  the  case  where  there  are  several  components  with  reac¬ 
tive  capacity  whose  constants  of  reaction  rate  exceed  that 
of  the  rest  of  the  components  by  a  considerable  margin,  the 
line  of  reasoning  is  analogous.  The  quasi-stationary  values 
in  this  case  are  determined  from  the  solution  of  a  set  of  al¬ 
gebraic  equations  obtained  by  setting  the  right-hand  parts  of 
the  corresponding  equations  in  set  (1)  to  zero.  The  times 
when  these  components  appear  in  the  quasi-stationary  condi¬ 
tion  is  inversely  proportional  to  the  corresponding  values  of 
Bk.  Since  the  constants  of  reaction  rate  incorporated  in 
br(t)  materially  depend  on  temperature,  the  time  of  appearance 
or  ck(t)  in  the  quasi-stationary  condition  under  various  con¬ 
ditions  of  the  course  of  the  reaction  may  change  considerably, 
decreasing  with  an  increase  in  temperature.  Prom  the  deter¬ 
mination  of  &  (t)  it  can  be  seen  that  the  magnitude  of  a,(t) 
influences  the  magnitude  of  corresponding  quasi-stationary 
value  ck(t). 

3.  The  algorithm  of  numerical  integration  of  the  set 
of  equations  of  chemical  kinetics  under  conditions  of  equa¬ 
tions  (4)  and  (  ™  ™* 

As  was  shown  earlier,  if  in  the  right-hand  side  of  the 
equation  for  any  concentration  there  are  one  or  several  mem¬ 
bers  containing  the  indicated  concentration  with  coefficients 
sufficiently  large  in  the  vicinity  of  the  modulus,  then  for 
a  short  time  interval  (compared  with  the  duration  of  the  reac¬ 
tion)  the  concentration  in  question  for  all  practical  purposes 
reaches  its  quasi-stationary  value  for  that  time.  It  follows 
from  equation  (3)  that  when  the  conditions  of  equations  (4) 
and  (5)  are  fulfilled,  the  concentrations  ck(t),  having 
reached  (to  the  degree  of  accuracy  required!  their  quase- 
s tat ionary  values  when  tk  =  (lnE\  /B^,  will  remain  in  their 
quasi-stationary  state  when  t  >  tk  also. 

However,  the  numerical  integration  of  this  set  of  non¬ 
linear  equations  which  is  under  examination,  including  in  it¬ 
self  one  or  several  equations  satisfying  the  conditions  of 
equations  (4)  and  (5),  leads  to  an  oscillation  of  the  solution 
around  its  quasi-stationary  value  for  practically  any  path 
of  integration  different  from  zero,  when  t  >  t  .  This  pheno¬ 
menon  is  connected  with  the  discreteness  of  counting  in  nu¬ 
merical  integration,  since  the  calculation  of  the  right-hand 
parts  in  equation  set  (1)  is  conducted  at  fixed  points  of  the 
integration  path.  For-  explanation  we  shall  introduce  the 
designation  zktt)  =  cj<(t)  -  H^(t),  Then,  subtracting  from 
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equation  set  (3)  the  equation  which  it  has  been  determined 
satisfies  the  concentration  ^(t) 

0 + MO. 

we  get,  using  the  Lagrange  theorem, 

tk  (<)  =  -  Mk  (<)  (<). 

where 

$6  (MO.  MO).  if  M0<M0= 
se(?» co,  MO).  if  mo<mo- 


When  the  conditions  of  equations  (4)  and  (5)  are  fulfilled, 
we  have  from  equation  (7)  that  a  small  increase  in  lz(t)| 
leads  to  a  large  change  in  ck(t),  as  a  result  of  which  )z(t)| 
decreases.  Hence  with  the  passage  of  time  z(t)  tends  toward 
zero.  However,  in  numerical  integration  with  an  interval 
.  zit^  l/B  (which  is  sufficiently  small  to  meet  the  condi¬ 
tions  of  equations  (4)  and  (5),  zJ t),  and  hence,  c,  (t)  change 
sign  within  the  limits  of  one  interval.  As  a  result  the  nu¬ 
merical  integration  of  equations  (2),  which  satisfies  the  con¬ 
ditions  of  equations  (4)  and  (5),  leads  to  an  oscillation  in 
the  solution  of  cMt)  relative  to  &  (t).  Decreasing  the  inter¬ 
val  of  the  integration  process  naturally  reduces  the  ampli¬ 
tude  of  the  oscillation.  However,  when  c,,(t)  appears  in  the 
quasi -stationary  state,  that  is,  when  z(t7  is  located  within 
the  limits  of  accuracy  for  the  calculations,  the  integration-- 
even  with  the  smallest  possible  interval — leads  to  oscillation 
whose  amplitude  increases  with  time.  Therefore,  the  numerical 
integration  of  this  set  of  equations,  one  or  several  of  whose 
equations  satisfies  the  conditions  of  equations  (4)  and  (5), 
leads  to  an. extremely  small  integration  interval  ( intentionally 
not  increasing  over  time)  and,  hence,  to  a  long  computation 
time  and  the  accumulation  of  error.  If  the  integration  is 
conducted  with  a  constant  integration  interval,  the  accuracy 
of  the  computation  quickly  falls,  and  the  oscillation  in  the 
solution  occurs  with  sharply  increasing  amplitude  leading  to 
computer  overload. 

As  is  known,  under  certain  circumstances  a  chemical  pro¬ 
cess  may  have  an  oscillatory  character  [source  However, 

if  it  were  possible  to  lower  the  amplitude  of  the  oscillation 
of  the  numerical  solution  caused  by  the  discreteness  of  the 
computation  o_n  the  computer,  the  amplitude  of  the  "real" 
oscillation  [source  83  caused  by  the  oscillatory  nature  of 
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chemical  process,  would  remain  constant.  Moreover,  the  amp¬ 
litude  of  the  oscillation  of  the  solution,  caused  by  the  dis¬ 
creteness  of  the  computation,  could  also  grow  without  limit, 
which  could  lead  to  violation  of  the  law  of  conservation  of 
matter  and  even  cause  the  appearance  of  negative  values  for 
the  concentrations.  Such  phenomena  do  not  occur  if  the  oscil¬ 
lations  of  the  solution  of  equation  set  (1)  are  "real.'1  Thus, 
the  different  behavior  of  the  two  types  of  oscillations  of 
the  solution  of  equation  set  (1)  in  principle  allows  us  to 
tell  these  oscillations  one  from  the  other. 

On  tne  basis  of  the  foregoing,  the  following  algorithm 
is  proposed  for  solving  problems  in  chemical  kinetics  when 
chemically  active  particles  are  participating  in  the  reac¬ 
tion.  The  integration  of  equation  set  (1)  in  this  case  must 
be  carried  out  in  two  stager. 

In  the  first,  transitional  stage,  where  the  concentration 
magnitude  materially  depends  on  the  selected  initial  condi¬ 
tions,  the  set  is  integrated  like  the  problem  of  Koshi  by  one 
of  the  methods  of  numerical  analysis.  The  duration  of  the 
first  stage,  as  was  pointed  out  earlier,  is  determined  by  the 
magnitude  of  the  corresponding  values  or  b^t).  Lin  practice, 
the  duration  is  of  the  order  of  10/min  B.,  if  there  are  several 
components  capable  of  reaction  whose  equations  satisfy  the 
conditions  of  equations  (4)  and  (5)J  Here  the  integration 
absolutely  must  be  carried  out  using  an  automatic  selection 
of  integration  interval,  but  the  initial  value  for  the  inter¬ 
val  must  be  taken  less  than  l/max  B^. 

When  it  is  taken  into  account  that  the  magnitudes  of  the 
concentrations  differ  from  each  other  by  many  orders  of  mag¬ 
nitude,  but  that  scaling  them  is  difficult  because  of  their 
rapid  change  over  time,  it  can  be  seen  that  the  integration 
must  be  conducted  to  a  given  relative  degree  of  accuracy;  that 
is,  to  a  certain  number  of  significant  figures  regardless  of 
the  magnitudes  themselves. 

The  final  goal  of  the  process  of  integrating  the  prob¬ 
lem  in  the  first  stage  is  the  arrival  of  the  components  capable 
of  reaction  into  the  quasi -stationary  state.  After  each  inte¬ 
gration  interval  a  comparison  is  made  between  c.  ( t )  and  &  (t). 
If  the  difference  between  the  solution  to  the  equation  having 
the  characteristic  examined  above  and  its  quasi -stationary 
value  for  any  moment  of  time — if  this  difference  is  within  , 
the  accuracy  limits,  then  we  can  go  into  the  second  stage  of 
the  integration.  Here  the  order  of  the  set  is  lowered,  and 
the  differential  equation  that  has  been  described  is  replaced 
by  an  algebraic  equation  (by  setting  its  right-hand  side  equal 
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to  zero).  In  the  general  case,  when  there  are  a  number  of 
elements  capable  of  reaction,  a  set  of  non-linear  algebraic 
equations  results  from  the  replacement.  As  will  be  shown 
later,  when  the  conditions  of  equations  (4)  and  (5)  are  met, 
such  a  replacement  does  not  introduce  divergence  and  does  not 
yield  an  accumulation  of  error,  but  instead  the  effective¬ 
ness  of  the  computation  is  sharply  increased,  enabling  us  to 
eliminate  the  difficulties  described  above. 

In  the  general  case  it  is  necessary  to  solve  the  set  of 
non-linear  algebraic  equations  obtained  in  the  second  stage 
by  integration  methods  at  each  referral  to  the  computer  pro¬ 
gram  for  the  right-hand  sides  of  the  equations  in  the  set  of 
differential  equations.  In  so  doing,  the  computation  time  • 
may  significantly  increase.  Hence  it  is  expedient  to  solve 
such  problems  on  high-speed  computers. 

4.  Evaluation  of  the  accuracy  of  the  algorithm. 

The  proposed  algorithm,  with  the  successive  replacement 
of  the  differential  equations  with  algebraic  equations  for 
the  several  components  with  high  reaction  capability  which 
satisfy  the  conditions  of  equations  (4)  and  (5)  when 
tk  =  Jln£|  /%»  insures  that  the  solution  of  equation  set  (1) 
is  obtained  to  the  desired  degree  of  accuracy. 

For  proof  of  this  we  shall  examine  the  change  in  deviation 
of  zk(t)  of  the  solution  of  equation  (3)  from  its  quasi-sta- 
tionary  state  when  t  >  t..  Taking  into  account  the  deter¬ 
mination  of  2fk(t),  the  equation  for  z,  (t)  from  equation  (3) 
may  be  written  in  the  form 

where,  in  accordance  with  our  earlier  examination, 

bvCtjC^).  It  follows  from  equation  sets  (1)  and  (2) 
that  when  expression  (4)  is  achieved,  \av(t)]  has  the  order 
a^b,, ,  and  It',,, |  the  order  b£.  Therefore < 

|**  +  A/***I<‘V*.  <*> 

where  is  some  value  which  is  the  upper  boundary  of  the 
expression  that  is  in  the  left-hand  side  of  equation  (3). 

If  the  condition  of  equation  (5)  holds,  ?,  . 

A  K 

We  shall  examine  the  integral  curves  (conation  S)  and 
the  straight  lines  z.  -  ±  ,  where  'l,  -  /\.  It  can  be 

seen  that  if  for  some  t*  >0  wo  have  z,.(t  (f  >*'H.  then 
z^'t*)  >  0,  and  the  integral  line  of  equation  (^J  may  only 


decrease  as  far  as  its  intersection  with  the  straight  line 
zk  =  R- .  If  z.  (t*)  >Rk,  then  z.  (t*)  >  0.  and  the  integral 
curve  of  equatxon  (8)  increases  until  it  intersects 
z  =  -R  But  if  lz  (t*))^.  Rw,  then  the  integral  curves  of 
equatioft  (8),  for  t^>t*,  can  not  intersect  the  corresponding 
straight  lines  2^  =  4  R^.  Henoe: 

{**  (f );  —  <  **  (*)  <  max  {tk  (ty,  fiky 

for  t  t*  and  for  t*  t^.  When  the  components  with  high 
reaction  capability  reachK(to  the  desired  accuracy  )  the 
quasi -stationary  values,  the  following  evaluation  occurs: 

min  {e;  /?,}  <  *,  (/)  <  max  {*;  fik}. 


Because  under  the  conditions  of  equation  (5)  Rk  is  quite 
small,  the  difference  between  the  quasi -stationary  solution 
and  the  real  solution  will  be  within  the  limits  of  the  accu¬ 
racy  of  the  computation,  when  t  >  t^.  Thus  the  proposed  algo¬ 
rithm  of  the  solution  of  a  set  of  differential  equations  de¬ 
scribing  the  kinetics  of  complex  chemical  reactions  when  com¬ 
ponents  of  high  reactive  capability  are  present  that  satisfy 
the  conditions  of  equations  (4)  and  (5)— this  algorithm  insures 
the  obtaining  of  a  solution  with  sufficient  accuracy  for  prac¬ 
tical  purposes. 

The  above  evaluation  of  the  magnitude  of  the  quasi-sta- 
tionary  concentrations ,  as  well  as  the  qualitative  evaluation 
of  the  accuracy  of  the  algorithm,  both  of  which  depend  on  the 
ratios  of  the  constants  of  reaction  rates  of  the  elementary 
stages  of  the  reaction  (that  is.  the  degree  to  which  the  con¬ 
ditions  of  equations  (4)  and  (5;  are  fulfilled) — these  evalu¬ 
ations  may  be  obtained  in  a  similar  way  in  every  concrete 
case.  It  is  evident  that  these  evaluations  are  excessive 
conservative  ,  Therefore  when  the  equations  of  balance  are 
done  (and  these  are  the  practical  standards  for  judging  the 
results  of  the  integration  of  equation  set  ( 1 ) )  it  will 
actually  be  done  with  greater  accuracy  than  the  evaluations 
obtained  above  would  indicate. 

We  note  in  conclusion,  that  in  the  case  of  reactions  of  the 
second  order  and  higher,  after  arrival  into  the  quasi-sta- 
t ionary  state,  it  is  essential  from  time  to  time  in  the  course 
of  the  solution  to  check  that  the  Conditions  of  equations  (4) 
and  (5)  are  being  met.  If  during  the  reaction  process  b. 
decreases  so  that  b^  <_  |ln£.l  /  At,  where  At  is  the  interval 
of  integration,  then  in  order  to  get  the  required  accuracy  it 
is  essential  to  replace  the  corresponding  algebraic  equation 
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with  the  initial  differential  equation,  etc. 


This  method  that  has  been  developed  by  us  above  will  be 
used  in  the  following  sections  of  this  chapter.  It  will  be 
applied  to  various  problems  of  the  kinetics  of  complex  chemi¬ 
cal  reactions. 

II,  A  calculation  of  the  kinetics  of  chlorination  of  ethy¬ 
lene  and  ah' "analysis  of  the  accuracy  o?  the  algorithm  of  the 
solution.  "  . 

The  complex,  multi-stage  reaction  of  the  chlorination  of 
ethylene  offers  considerable  difficulties  for  investigation. 

In  this  section  we  obtain  the  functions  c*  =  c.(t)  for  con¬ 
centrations  of  intermediate  and  final  proaucts^of  this  reac¬ 
tion  by  the  mechanism  of  combination  in  the  temperature  inter¬ 
val  600  -  1200°K  and  for  differing  initial  ratios  of  CI2  :  CgH, 
An  appropriate  analysis  of  the  accuracy  is  conducted  for  the  ^ 
results  of  a  numerical  integration  of  a  set  of  differential 
equations  of  the  kinetics  of  the  chlorination  of  ethylene. 

This  analysis  also  has  a  more  general  significance  for  the 
nature  of  the  method  of  numerical  integration  of  kinetic  equa¬ 
tions  of  complex  chemical  reactions  that  was  developed  in  sec¬ 
tion  T-3  of  this  paper. 

On  the  basis  of  the  data  now  available  it  is  considered 
that  substitutional  chlorination  may  proceed  in  two  direc¬ 
tions  (sources  8,9]  : 

<*H*  +  Cl,  C,H*C1„ 

0,11*01,  -V  C,H,C1  +  HC1 ,  (9) 

C,H«  +  Cl,  —  C,H,C1  -f  HC1.  (10) 

In  the  first  case  vinylchloride  is  formed  from  the  pro¬ 
ducts  of  combination  (C2H^C12),  and  in  the  second  case  there 
is  a  direct  replacement  process  [source  10]  .  On  the  basis 
of  present  computations  the  first  of  these  mechanisms,  which 
is  considered  the  basic  mechanism  by  many  authors,  is  assumed. 

A  reaction  process  v/as  constructed  for  mechanism  (9)  that 
was  considered  to  be  one  of  the  possible  models  for  the  chlo* 
rination  of  ethylene. 

1.  Formulation  of  the  -problem. 

V/e  have  examined  the  following  reaction  process,  in 
which  elementary  stages  known  in  the  literature  are  used 
[sources  11-14]: 


1)  Cl  4-C,H*-^CaH4CI 

2)  C,H4C1  -»a  +  C,H4 
3}  0,^01  + Cl -i-CjHiC], 

4)  C*H*C1  +  HC1  —  C,HSC1  +  Cl 

5)  C.H.C1  +  Cl  —  C,HtCl  f  HCl 

6)  C,H4C1  -f  Cl,  -  C,H4C1,  +  Cl 

7)  c,h4ci, + a  -*  c,h4ci + ci, 

8)  q,II«Cl,  +  Cl  ->C,!IaCl,  +  HC! 

9)  C,HaCl,  -f  HCl  —  C,H4C1,  +  Cl 

10)  CjHjCl,  +  Cl,  -»  C,HSC13  +  Cl 

11)  C,H0C11  +  C1->C,HSCI,  +  C1, 

12)  C,HaCl,  +  Cl-*C,H3ClJ 

13)  C,HjCi,  -►  C,HaCl  -f-  Cl 

14)  C,H,C1  +  C1— C,HaCl, 

15)  C1,  +  M-*CI  +  C1  +  M 

16)  a  +  Cl+M— Cl,  +  M 

17)  C,HtCl  +  C;H4CI  -»npo;iyKTH 

18)  Q,H,C1S  -f  CI-+  C,H,C1S  +  HCl 


(It) 


*?he  authors  thank  V.  Ya.  Stern  and  A.  P.  Revzin  for  help 
given- in  the  construction  of  this  reaction  process. 


Since  further  chlorination  does  not  interest  us,  we  limited 
ourselves  to  the  reactions  listed  above. 

Corresponding  to  the  reaction  process  (11)  was  constructed 
a  set  of  12  non-linear,  ordinary  differential  equations  (12) 
which  were  numerically  integrated  on  the  "Strela"  (_' Arrow '1 
computer  of  the  computer  center  of  Moscow  State  University. 

^hs  integration  was  carried  out  using  the  algorithm  of  inte¬ 
gration  of  a  set  of  equations  of  kinetics  of  complex  chemical 
reactions  which  was  examined  earlier  in  this  paper  in  section 
1-3. 


-  16  - 


Wi  +  Wi  —  kyfiiCx  +  kac&  —  *ii<Vo4- 

jf  —  —  4-  fat  —  kte4,et  +  k^et  —  kxfijCt  +  htfitPi  — 

—  *7<V» — V*e»  +  *«e»e*  4  *io<Vi  —  Aruc^,  —  A-Wc»c,  + 

4*  kufix  —  kuevfi*  4*  24i*CiCo  —  2fcucjcg  —  iuV> 

*j|j  ~ — AjCjC, 

^  =  V»c3  —  *«e«  —  ^*c*ei  —  V«c«  4  ArfjCt  —  V|fi  4 

+  **»««  — 

37  “  V«c,  4  Va  —  fc?e*c*  —  V*c«  4  W* 

—  fr«c4c«  +  ^CjC,  +  *4c4e,  —  Aw,  +  A-jgCjc* 


77  =  ktctct —  he*, 

“  feVi  ^§*1  kia£$Ci  4*  ^uc«f»  —  ^iACi — ii}C|  4* 

4*  kuCj/fix 

*■*  *ttC*ex — AuCgCj  4-  ^ijCiCj  —  *u<y» 

7T  =  *jC$  —  A*«CloCl 

IT  -  *«cl 

“jjT  “ 


>(12) 


In  set  (12)  the  following  designations  are  hereby  adopted 
for  convenience: 


)  Mopene 

KQJyl 

(^Beaeotao 

^obhcbt- 

pwxmj 

<3. 

CgHjCl 

•t 

a 

C»HjCl* 

•• 

C.B. 

H 

CiBjCIj 

«• 

c,H*a 

€4 

C,H>C1 

CjHiCli 

<• 

©npoflyimr 

MHHiUtlaTOi 

PMKOJH 

Cu 

Ha 

«• 

c«H«a> 

*U 

Key: 

1,  Substance 

2.  Concentration 


3.  Substance 

4.  Concentration 

5.  Products  of  the  17th  reaction. 
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In  the  case  where  this  process  occurs  in  a  current,  equations 
of  hydrodynamics  should  be  added  to  the  set  of  differential 
equations  of  kinetics  (12)  and  the  set  thus  obtained  should 
be  solved  concurrently,  as  was  done  for  instance  in  source 

for  the  process  of  the  conversion  of  methane  to  acetylene,- 

For  the  initial  conditions  of  the  equation  set  (12)  we  use 
the  initial  concentrations  of  ethylene,  molecular  chlorine, 
and  elemental  chlorine,  where  the  initial  value  for  the  con¬ 
centration  of  elemental  chlorine  is  found  by  means  of  a 
balanced  degree  of  dissociation  of  molecular  chlorine  at 
the  moment  t  =  0. 

The  values  of  the  constants  of  reaction  rate,  which  are 
necessary  for  numerical  integration,  were  computed  (with  the 
exception  of  the  constant  for  the  dissociation  and  reassocia¬ 
tion  of  chlorine)  by  means  of  the  well-known  formula: 

{—E{  l  RT),  (13) 

where  the  factor  before  the  exponential,  kg.,  and  the  energy 
of  activation,  E.  were  taken  from  sources  (J2,  16,  ly  which 
have  the  data  which  is  most  complete  and  mutually  substan¬ 
tiating.  These  data  are  listed  in  Table  1,  Because  data  are 
lacking  for  the  higher  temperatures,  the  values  indicated  for 
the  constant  of  reaction  rates  for  the  1200°K  region  were 
extrapolated  by  us. 

Table  1. 

Values  of  the  energy  of  activation  and  pre-exponential 
factors  for  the  elementary  stages  of  the  reaction  of 


< 

tffi 

i 

,<3> 

- ST" 

Ei 

*oi.\iw3 

&{  ,  |*C  K«V-HOU«fc  \ 

1 

0 

10 

6,3-10% 

1 

2 

7,9-iO11 

23,6 

11 

2-10*% 

20,6 

3 

2-10'% 

0 

12 

2-10*% 

0 

4 

10% 

9.1 

13 

6,3«10** 

23,8 

5 

3,2-10'% 

1.5 

14 

1,6.10'% 

0 

6 

2,5-10% 

1 

17 

1,3.10*% 

0 

7 

2*10*% 

21,3 

18 

3,2.10*% 

3,5 

8 

2,5.10*% 

3 

6 

7,9-10% 

10,2 

Key: 

1.  l/sec.  3»  l/sec. 

2.  kilocal. /mole  4.  kilocal./mole 

The  reaction  process  (11)  also  embodies  reactions  of 
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dissociation  and  recombination  of  chlorine. 

A  number  of  authors  have  shown  in  experiments  using  impact 
tubes  (sources  18,  19]  in  the  temperature  interval  from  1500 
to  2500°K  that  the  experimental  data  on  the  constants  for  the 
dissasBOOlation  of  chlorine,  for  corresponding  choices  of 
n,  kQ,  and  E,  are  equally  well  described  by  either  of  the  two 
following  formulas: 

=  a  YT  (Da/RTf  oxp  (~D6/RT),  (14) 

kd>\vT  *° exp  E/RT)'  (15) 

where  n  is  the  numbered  parameter.  Usually,  n  is  taken  to  be 
n  «•  2  "^sources  18,  19]  . 

As  far  as  the  numerical  value  of  the  energy  of  activation 
is  concerned,  different  authors  give  values  for  this  temperature 
interval  that  do  not  agree  with  each  other— from  40  to  48  kilo- 
cal./mole.  [sources  17  -  19]. 

For  the  integration  of  equation  set  (12)  in  this  study, 
the  Constants  of  disassociation  and  recombination  were  com¬ 
puted  according  to  the  following  formulas  [source  18]  : 

klb  =  1,1.10*  T'/‘  (D0/RT)*  exp  {-(Da/RT)],  (16) 

ku  =  3,6-10'  (D0/RT)*/[i  -  exp  (-813/f)].  (17) 

2.  The  algorithm  of  bhe  solution  of  equation  set  (12). 

From  what  has  been  written  above  it  follows  that  the  con¬ 
stants  of  rates  of  reaction  change  within  wide  limits  over  dif¬ 
ferent  stages  of  the  reactioiiounder  investigation  (for  example, 
at  temperature  600°K  from  10"  to  10°).  As  was  pointed  out 
earlier  (middle  of  section  T-2),  this  wide  range  of  values 
creates  considerable  difficulties  for  the  numerical  integra¬ 
tion  to  a  desired  degree  of  accuracy. 

Tn  the  case  of  equation  set  (12)  the  components  whose 
equations  satisfy  the  conditions  of  equations  (4)  and  (5),  are 
C2H/jCl,  028*3012,  and-  also  under  certain  conditions  Cl  (our 
designations  c^,  eg,  and  C2  correspondingly) . 

In  the  temperature  range  from  600  to  1200°K  the  concen¬ 
trations  c/+(t)  and  eg ( t)  reach  their  quasi -stationary  states 
with  values  c!<:(t)  =.  ai^tJ/bjd t)  (k  =  4,3)  correspondingly ,  not 
later  than  from  1CT  to  lO"'8  sec.  after  the  start  of  the  reac¬ 
tion,  and  in  agreement  to  what  was  written  above,  (section  1-2), 
the  concentrations  shall  intentionally  be  allowed  to  stay  in 
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the  quasi-stationary  state. *  It  is  here  evident  from  equation 
set  (12)  that  c1(0)  =  0,  c2(0)  <  0,  c^(0)<.0,  6^(0)  >0. 

*As  was  shown,  under  certain  conditions  this  is  also  justified 
for  Cl. 


The  integration  of  equation  set  (12),  as  was  shown  earlier 
(section  1-3),  must  be  carried  out  in  two  stages.  The  initial 
interval  of  integration  must  be  intentionally  not  more  than 

mini*  (0)  (*  =  2, 4. 8), 

that  is: 

“*nfeJ(or'  i'  £}• 

In  this  case  c^(t)  and  cg(t)  enter  the  quasi-stationary 
state  almost  simultaneously,  Because  bjJt),  when  k  =  4,8,  (the 
lower  evaluations  of  which  are  correspondingly  equal  to  k  and 
k.-J  are  closed  to  each  other.  Since  c  (t)  and  (L(t).  ^accor¬ 
dance  with  equation  (12),  do  not  clearly  depend  on  c^Ct)  and 
Co ( t)  respectively,  in  this  case  5>,(t)  and  c_(t)  may^be  found 
independently  of  each  other  from  tne  equations: 

cg(t)  =  0 
and  £“(t)  =  0 

Besides  this,  the  condition  c.  (t)  =  c.(t)  =  0  allows  the  equa¬ 
tion  for  co(t)  to  be  brought  into  a  more  convenient  form  for 
numerical  integration: 

c%  (t)  =  2  (kiiPj c$  —  k AjjCoC,)  —  ki$CfCt, 


where  the  magnitude  of  the  coefficient  for  c?  decreases  signi¬ 
ficantly  compared  to  equation  set  (12).  ^ 

A  check  showed  that  when  equation  (4)  was  solved  using  the 
proposed  algorithm,  the  equation  of  balance  is  fulfilled  to  the 
desired  accuracy,  which  is  determined  by  the  relative  degree  of 
error  inherent  in  the  method  of  numerical  integration  over  the 
whole  integration  process.  The  computation  time  machine  time 
required  for  one  variant  of  the  Hunge-Kutt  method,  using  an 
automatic  chpice  of  interval  of  integration  and  a  relative 
error  of  10"  and  10-5 f  is  about  one  hour  on  a  "Strela-^1* 
(ArrowJQ  computer  with  a  computer  speed  of  about  2000  opera¬ 
tions  per  second. 
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3.  An  evaluation  of  the  accuracy  of  the  results  of  the 
numerical'  integration  of  equation  set  (12) 

We  shall  show  that  in  the  temperature  range  being  examined 
(600  -  1200°K)  the  algorithm  of  the  solution  of  equation  set 
(12)  insures  an  absolute  accuracy  not  lower  than  plus  or  minus 
lO-o. 


Beforehand  we  shall  obtain  the  upper  evaluation  of  the  mag¬ 
nitude  of  the  concentration  of  atomic  chlorine  as  a  function 
of  time  in  the  given  temperature  r^nge.  For  this  purpose  we 
shall  examine  the  equation  for  Cl  in  the  form 

i,  (/)  =  /  (clt  eu)  —  et  tk&  +  k,ej  +  (k,  +  k$)  c,  + 

+  ktfr  +  kuelB  +  2  4-  *i c,  +  (&u  +  *u)  c»l  +  c,  (kt  + 

4-  kfCg  4'  ktfj)  4-  c%  (ku  4"  kjCt  4"  AjoCi)  .4*  2fcj»eic#.  (18) 

We  shall  find  a  cjj?  such  that  for  c2  and  for  any 

other  possible  values  of  the  other  concentrations,  the  follow¬ 
ing  holds: 

1  (*JI  *»>•••!  ^||)  0. 


This  means  that  the  reaction  being  examined  may  proceed  with 
an  increase  in  the  concentration  of  atomic  chlorine  only  if 
c2  <  c£. 


As  follows  from  equation  (18),  the  positive  terms  do  not 
clearly  depend  on  c.,  c?,  eg,  and  c1Q,  We  shall  hold  the  values 
of  the  other  concentratfons'constant  and  write  the  equation  of 
balance  for  Cl  in  the  form: 


2e«  4-  4-  3e»  4-  Cj«  =  a  ~  2fj  —  c,  —c,  —  c,~  2c,  — 2cu— 3 clt, 

(19) 


where 


a  =  2  ct  (0)  +  et  (0). 

It  is  evident  that  the  function  f  will  reach  a  maximum  if,  for 
arbitrarily  fixed  values  of  the  indicated  concentrations,  the 
negative  terms  in  equation  (18)  are  at  their  minima.  This  con¬ 
dition  is  fulfilled  if  the  expression  obtained  from  equation 
(18)  and  presented,  taking  equation  (19)  into  account,  in  the 
form 

=  k,c,  4-  -g-  (k-i  4-  k,)2c,  +  kucl9  4-  --  {ku  4-  ku)  3c», 

is  at  its  minimum.  In  view  of  the  fact  that,  in  accordance 
with  Table  1,  the  following  inequality  holds: 
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•j-  (^'i*  +  ^11)  'C  m'n  ~2  (ki  4‘  A>),  Am|  , 

El  reaches  its  minimum,  (as  a  consequence  of  equation  (19), 
if 

Cio  -  0.  *  (20) 

Thus  the  condition  of  equation  (20)  is  essential  if  func¬ 
tion  f  is  to  reach  its  maximum  value.  Then  expression  (18) 
takes  the  form  "i 

/  M=  —  et  +  (An  +  ^u)  c»  +  l,c,)  -f 

+  c*  (*t  +  Ket  +  foi)  +  c*(kit  +  +  2  A^gC,.  (2i) 

Further  search  for  the  majorant  of  c2(t)  wili  be  accom¬ 
plished  by  means  of  varying  c.,,  c^,  eg,  of.,  and  c12,  which 
heretofore  were  held  constant! 

Let  us  look  at  the  contents  making  up  c.  and  c ,  in  for¬ 
mula  (21),  taking  into  account  expression  (26),  and°for  the 
values  of  Cn,  c^,  c  ,  c11f  c12,  and  being  held  arbitrarily 
constant.  Then  frofl  the  balance  in  equation  (19)  we  have 

2c,  +  e«  *■  a  —  3c»  —  eK  —  2c*  —  2eu  —  3cM  —  et.  (22) 

Consequently  the  expression  entering  the  right-hand  side  of 
equation  (21)  and  written  taking  equation  (22)  into  account, 
in  the  form 

fttct  d"  kfC$)  -f-  c,  (ktci  4"  kjgCi  *j-  2A^,c0), 

reaches  its  maximum  when  ct  =  0,  since  k<  >>  2k^,  2k_, 

and  also,  in  accordance  with  equation  (15),  k^ c 0.  Conse¬ 
quently  J 

■j  (A* 4  +  +  2^)  >  ktc,  4-  Ay,. 

From  the  above  it  follows  that  when  expression  (20)  is  satis¬ 
fied,  (regardless  of  the  values  of  the  other  concentrations), 
the  condition 

e*  "  0  (23) 

leads  to  the  maximum  value  of  c2(t)  for  the  given  state  of  the 
system.  As  a  result  of  formula^ ( 21 )  we  get 

(0  ^  ct  (^i  4~  kffii  kfC^j  -)■  c,  (Ajj  4"  —  k^c^  4* 

4-  2inCo^i  —  2kitC<£i  —  (Ajg  4"  kn)  CtCt  —  (24) 
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holds : 


It  is  evident  that  when  CgCO),  the  following 


(0  —  klt  4  (0  ^  o>  t  ^>0. 


Thus,  when  c?  =  c|,  the  only  terms  in  the  right-hand  side 
of  formula  {24)  which  may  be  positive  are  R_  =  kg  +  k^c,  -  k-c?i 
and  Rh  =  +  kinc1  -  k10c?).  In  accordance  With  Table  1; 

k?  >  k^,  k6>  k1Q,  and  kKfc  k.2  independent  of  the  teroperatur 
Hence,  If  Ro  ^  O^then  Rz,  <  0  aM  f(c  ,  c2,  „)  \  0. 

Consequently  a  maximum  of  X  i-3  possible  when^R-,  /  0.  Then 
R- >  Rn  and  c^R-  +  ZcaRn  will  not  be  large,  if^for  the  arbi¬ 
trary  values  ofKhe  other  concentrations  (taking  into  account 
equations  19,20,  and  23),  we  take 


Thus  from  expression  (24),  taking  into  account  equation  (26), 
we  get 

/  (^i i  tit  #*»•••!  ^it)  ^  (^*j  4  ktex  -  kfCi)  4 

+  2  ~«fri«eal)  e«  —  (^i*  +  ktl)  etfit  —  kx cft.  (27) 

Throwing  out  the  last  terras  of  equation  (2?)  (here  the  right- 
hand  side  of  equation  (27).  in  agreement  with  equation  (25), 
is  not  allowed  to  decrease)  we  finally  have: 


/  (®1*  .  •  i  £ll)  ^  0 


with 


4  =  max  (c,  (0),  c,}, 


where 

r  *•  +  *«ej 

When  there  is  a  surplus  of  ethylene,  c*L  =1,  and  since 
c.  (t)  <L  c. ,  (OK  <*■  /2 ,  c?  5*10“-;!  for  the  temperature  in¬ 
terval  from  600  to  1000°K7  Inasmuch  as  the  degree  of~disasso- 
ciation  in  this  interval  changes  from  10“  to  2.5'10~'> , 

Cg  y  anc^  the  evaluation  we  are  seeking  has  the  form 

«»(0)<c»(0  <5.10->.  (28) 

The  evaluation  obtained  in  equation  (28)  is  sufficiently 
effective  for  the  high  temperatures  within  the  limits  of  the 
temperature  interval  under  examination.  But  for  a  lowered 
temperature  the  evaluation  of  equation  (28)  becomes  substan- 
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tially  increased;  therefore  for  the  temperature  range  from 
600  to  800°  K,  (degree  of  dissassociation  of  which  is  10“B— 
10“5)  it  is  expedient  to  use  the  following  evaluation.  Prom 
formula  (12)  we  have 

(0  +  ( t )  +  c,  (t)  = 

**  —2  +  kljc\  -f-  kuSt et  +  ^  (29) 

tfe  shall  designate  c?(t)  +■  c.  (fc)  c2(t)  =  4>  (t),  whereby, 
since  0^,(0)  =  cg(0)  =  0,  we  nave  4>  (8)  =  c?(0).  Then  from 
equation  (29)  it  follows  that  <s>(t)  4  2k.c.';CjC0  for  any  t>  0. 
Hence  we  have 

T  (i)  <  2kltCle0t  +  (|  (0). 


Since  for  temperatures  up  to  600°K  the  rate  of  disassociation 
kir<,  10“?,  and  the  balanced  values  of  c2(0)  ^10“°,  for  a 
time  of  reaction  of  the  order,  of  one  second  we  have,  finally, 

(for  ^  ^  1,  cQ  -C  0.05); 

8  (0  +  *4  «  +  <?.  (0  <  2*„c*  +  C.  (0)  *  10'*. 

Taking  into  account  that  the  values  for  the  concentrations 
can  not  be  negative,  it  follows  that 

c*  (0  <  10-  (*  =  2,4,8).  (30)' 

At  800°K  c^Ct)  10“6.  It  is  evident  that  for  high  tem¬ 
peratures  this  evaluation  becomes  rougher  than  that  of  equa¬ 
tion  (28). 

The  rather  rough  upper  evaluation  of  the  magnitude  of 
c?(t),  carried  out  above,  permits  the  evaluation  of  the  magni¬ 
tude  ck(t)  (k  =  4,8)  in  the  reaction  process.  From  formula  (4) 
it  follows  that  ak(t)  (k  =  4,8)  is  proportional  to  c2.  In 
accordance  with  Table  1,  in  the  temperature  range  from  600°  to 
1200°K,  the  constants  of  reaction  rates  going  into  ak(t)  do  not 
exceed  109  in  the  dimensionless  form.  Consequently,  for 
temperatures  from  600  to  1200°K,  ak(t)  does  not  exceed  10  and 
5*10“  Psicq  respectively.  In  this  range  and  for  these  condi¬ 
tions, ^^Tt)  (k  =  4,8)  is  not  lower  than  2,5'10^  and  4*10y, 
respectively. 

As  a  result  in  the  temperature  range  being  examined  the 
upper  evaluation  for  the  magnitude  of  the  quasi-stationary  values 
of  the  chemically  active  elements  in  the  reaction  process  has  the 
following  form: 

«»(0<10-*,  <>0.  *  =  4,8.  .  (31) 
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It  is  evident  that  for  temperatures  of  less  than  1200°K  the 
evaluation  of  equation  (31)  becomes  quite  increased. 

With  the  help  of  the  results  that  have  been  obtained  we 
shall  find  a  qualitative  evaluation  of  the  accuracy  of  the  pro¬ 
posed  algorithm  of  the  solution  of  equation  set  (12),  this  so¬ 
lution  being  connected  with  the  replacement  of  the  differential 
equations  by  algebraic  equations  for  the  components  capable  of 
reaction  (Cl,  C2H/1CI  and  C2H2C12)  for  all  t  >  10"5(in£.\  . 

Here  the  same  designations  will  be  used  that  were  usea  in  sec¬ 
tion  11,1,  where  a  proof  was  offered  for  the  accuracy  evaluation 
in  its  general  form. 

We  shall  evaluate  the  magnitude  N,  (t)  (see  formula  8, 
section  I).  Prom  equation  set  (12)  itKcan  be  seen  that  the 
only  term  of  the  first  order  not  containing  c.(t)  (k  »  2,4,8) 
is  ki£cic0.  This  term,  which  characterizes  tne  disassociation 
of  chlorine,  does  not  exceed  2. 5»  10"-^  in  the  temperature  range 
from  600  to  1200°K  and  is  quite  small  compared  with  BSJ'  10^. 

Only  concentrations  of  the  reactive  components  c^,  and  eg  go 
into  the  other  first-order  terms.  In  addition,  In  the  second- 
order  terms,  at  least  one  of  the  concentrations  coincides  with 
0304  or  eg.  Thus,  in  agreement  with  the  evaluations  of  equation 
(2P;  and  equation  (31),  c(t),  j  =  1,2,.*, 12,  are  at  least  tforee 
orders  of  magnitude  less  than  the  maximum  for  the  constants 
of  reaction  rates  in  the  positive  terms  in  the  right-hand  sides 
of  the  corresponding  equations.  Here  the  derivatives  of  the 
concentrations  of  the  reactive  components  after  entry  into  the 
quasi -stationary  state  are  quite  small.  This  circumstance  enables 
us  to  obtain  an  evaluation  for. the  magnitude  of  a.  (t),  &Jt) 

(k  =  4,8)  and  also  for  a2(t),  b2(t)  when. there  isKa  surplus  of 
ethylene  present  and  t  >  t  .  Actually,  a^(t)  (where  k  =»  2,4,8) 
is  composed  of  the  final  sum  of  terras  having  the  form  k*,  Ci 
and  k^  (cmcQ  +  ,  which  are  obtained  by  differentiating*^ 

the  respective  first  and  second  order  terms.  It  is  evident 
that  the  first  of  the  indicated  terms  are  small,  since  the 
corresponding  components  are  capable  of  reaction. 


Let  us  assume  that  c  has  little  capability  of  reaction 
(then  c  is  quite  reactive).  Hence  in  the  worst  case  from  the 
standpoint #of  evaluation  (when  T  =  1200°K)  we  shall  obtain 
lUc  c  t  c  Cjj)  10”°k.kr,  where  kr  is  the  greatest  of  the 
constants  or  rates  of  reaction  in  the  right-hand  side  of  cn(t). 
If  cn  and  cra  are  capable  of  reaction,  then  the  corresponding 
term  will  be  quite  small.  As  a  result  we  have 


«»  to  <  10**. 


> 
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In  an  analogous  manner,  since  b.  (t)  is  composed  of  terms  of 
the  form  k^cLft),  where  c  may  nave  a  small  capability  of 
reaction,  we^have  p 

Hence  in  the  worst  case  we  -obtain 

since  for  a  temperature  of  1200° 

bk  ( t )  >  10*. 

For  t>  tQ  the  error  in  replacing  the  concentrations  of  the 
indicated  components  with  their  quasi-stat ionary  values  satis- 
fies  the  condition 


xk  (I)  <  max  {«, 


(32) 


where  Rk  <  10“^  at  1200°K.  As  follows  from  equation  (30), 
for  lower  temperatures  is  still  smaller. 

Thus,  the  proposed  algorithm  describes  the  kinetics  of 
chlorination  of  ethylene  and  insures  the  obtaining  of  a  solu¬ 
tion  with  an  accuracy  sufficient  for  practical  purposes, 

4.  A  discussion  of  the  results  of  the  numerical  inte¬ 
gration  of  equation  set  (12). 

Numerical  integration  of  equation  set  (12)  was  carried 
out  for  various  initial  conditions  of  the  ratio  of  the  concen¬ 
tration  of  chlorine  to  the  concentration  of  ethylene  in  the 
temperature  range  from  6 00  to  1200°K,  The  corresponding  mathe 
matical  experiments  were  carried  out  on  a  computer  for  evalu¬ 
ating  the  influence  on  the  solution  of  equation  set  (12)  of 
the  following:  changes  in  the  constants  of  reaction  rates 
kja  and  changes  in  the  reaction  processes  (equation  set 

(11))  by  means  of  excluding  individual  stages  (particularly, 
reactions  17  and  18  of  equation  set  (11)), 

Figures  la  and  lb  show  in  semi-logarithmic  scale  the 
results  of  the  numerical  integration  of  equation  set  (12) 
without  taking  into  account  reactions  17  and  18  in  reaction 
process  (11),  and  with  the  full  set  of  12  differential  equa¬ 
tions  [.equation  set  (12)3,  respectively.  The  calculations 
carried  out  showed  that  the  exclusion  of  reactions  1?  and  18 
from  the  reaction  process  (11)  had  insignificant  influence  on 
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the  time  dependence  of  the  concentrations  of  the  majority  of 
the  components  taking  part  in  the  reaction.  What  influence 
there  was  exerted  to  a  greater  degree  on  the  concentrations 
of  the  components  which  are  directly  connected  with  reactions 
17  and  18. 


Key: 

1.  c.  in  volumetric 
snares 

2.  c  in  volumetric 
shares 

3.  t  in  sec. 

4.  £  in  sec. 


Figures  la,  lb.  The  dependence  on  time  of  the  concentrations 
c<  of  various  components  of  the  reaction  a)  excluding,  and 
b;  including  reactions  number  17  and  18  of  equation  set  (11) 
in  the  reaction  process. 


Initial  ratio 

1)  lOClpJ  l  3)  10lC2Hj7  r5)  ioibXci2i ; 
9)  CgHoCl^  ;  10)  lOEXci}.  in  Figures 


to  Gig  to  Cr 


H,  =  0.3 


1 — 3  an 


0.7;4T  =  120QjK; 

8) 


a  I?- 


.  .  5)  10  ii;2“4VJ-2-i  *  *v  i«2li3'' 

^HnCl}..  in  Figures  la,  lb  the  numbers 
9  on  the  curves  correspond  to  the  numbers  of  the 


components  in  the  table  immediately  following  equation  set  (12) 


In  order  to  study  the  influence  of  the  [initial]  conditions 
of  the  proportions  of  the  concentrations  on  the  kinetics  of 
the  chlorination  of  ethylene  a  numerical  integration  of  equa¬ 
tion  set  (12)  was  carried  out  for  various  ratios  of  chlorine 
to  ethylene.  (3ee  Figure  2) 
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Key: 

1.  c.  in  volumetric 
snares 

2.  t  in  sec. 


Figure  2.  The  dependence 
on  time  of  concentrations 
of  chlorine,  ethylene,  and 
vinylchloriae  for  a  tempera¬ 
ture  of  1200°K  and  various 
initial  proportions  of  the 
c  oncentrat i ons . 


C2%  - 


Cl 
Cl 

ci  :  c;h£  ■ 

1)  Cl2f  v3) 


C2H4  " 


0.3  : 

0.5  : 
0.?  : 
C2H4; 


0.7  (solid  lines) 

0.5  (mixed  dotted  and  dashed  lines) 
0.3  (dotted  lines) 

10)  c2h3ci 


V/e  examined  the  case  where  there  was  a  surplus  of  chlo¬ 
rine-"’  present,  the  case  of  a  surplus  of  ethylene,  and  also  the 
intermediate  case  between  these  two  cases.  In  practice,  the 


*0f  course,  such  concentrations  of  Cl_  are  examined  only  for 
the  purposes  of  explaining  the  basic  tendencies  of  the  pro¬ 
cess;  in  practice  a  phenomenon  having  an  explosive  nature 
occurs  for  such  concentrations. 


chlorination  of  ethylene  is  usually  carried  out  with  an  initial 
surplus  of  ethylene  in  order  to  maintain  isothermal  conditions 
and  to  avoid  having  to  purge  the  reaction  products  of  free 
chlorine. 

As  can  be  seen  in  Figure  2,  in  the  case  where 
Zl-2  :  'W"'h  -  0.3  •’  0.7,  for  a  time  span  of  about  one  second, 
the  chlorine  is  practically  expended;  thereupon  the  c  ■>ncentra- 
tions  of  vinylchloride  and  ethylene  do  not  change  further.  It 
’ s  evident  that  the  graph  lines  for  concentrations  of  ethylene, 
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vinylchloride,  and  a  number  of  other  components  become  asymp¬ 
totic  because  the  chlorine  has  been  expended.  In  order  to 
exclude  this  factor,  a  variant  where  Cl2  :  CLIfy  =  0*7  :  0.3  ■ 
was  calculated.  As  can  be  seen  from  Figure  2,  the  kinetics 
of  the  process  in  this  case  differ  markedly  from  that  of  the 
former  case:  the  concentration  of  vinylchloride  reaches  a 
maximum  value,  after  which  its  intensive  decomposition  occurs. 
At  the  moment  when  the  concentration  of  vinylchloride  reaches 
its  maximum  value,  the  concentration  of  ethylene  is  about 
0.06  and  that  of  chlorine  is  about  0.40  (concentrations  of  the 
substances  are  given  in  volumetric  shares). 

Thus  in  the  case  where  Cl2  :  CpH,  =  *  °*3>  as  could- 

have  been  expected,  the  process  tends4 toward  a  more  intensive 
chlorination  than  in  the  case  Cl2  :  C  H,  =  0.3  :  0.7;  the 
appearance  of  vinylchloride  is  decreased,  and  the  concentrar 
tion  of  HC1  continues  to  increase  even  after  the  concentration 
of  vinylchloride  has  reached  its  maximum.  Also  calculated  was 
the  intermediate  case  where  Cl2  :  C2%  =  0.5  :  0.5,  As  can  be 
seen  in  Figure  2  the  concentration  of  vinylchloride  reaches  a 
maximum  value  of  about  0.25  and  thereafter  practically  does  not 
decompose  at  all. 

Also  investigated  was  the  question  of  the  influence  of 
the  constant  of  reaction  rate  of  reaction  number  13  in  the 
reaction  process  described  by  equation  set  (11)  ,  which  is 
the  only  reaction  in  the  reaction  process  we  have  chosen  that 
leads  to  the  formation  of  vinylchloride.  To  this  purpose 
equation  set  (12)  was  numerically  integrated  for  a  temperature 
T  =  600-K,  an  initial  ratio  of  Cl2  :  C?Hj,  =1  :  10,  and  for 
values  of  Efjo  10  and  100  times  larger  than  the  values’  of  ki-> 
corresponding  to  600 °K.  ' 

Figure  3  shows  the  appearance  of  vinylchloride  at  time 
t  =  0.3  sec.  as  a  function  of  the  value  of  k  /k-,.  It  is 
easy  to  see  from  this  figure  that  the  furtheJ*3 increase  in 
value  of  this  constant  does  not  lead  to  a  substantial  increase 
in  the  appearance  of  vinylchloride. 

Figure  3.  The  dependence  of  the 
appearance  ofj/inyichloride  on 
the  value  of  K  /k^, 

t  -  600°K;  Cl  :  CH  -  1  ;  10 

Key;  1)  c£J  in  volumetric 

^  shares. 


[CjHjCi],  0618MH.  ftvwfl) 


In  connection  with  the  great  sensitivity  of  the  process 
to  changes  in  the  constant  of  reaction  rate  of  disassociation, 
the  question  of  the  influence  of  changes  in  magnitude  of  this 
constant  was  studied.  There  are  no  exact  measurements  of  this 
constant;  moreover,  the  possibility  cf  the  walls  having  a  sig¬ 
nificant  influence  on  its  magnitude  must  be  taken  into  account. 
Therefore  calculations  were  carried  out  for  various  values  of 
the  energy  of  activation  of  disassociation  and  various  values 
of  the  constant  of  reaction  rate  of  disassociation. 


Figure  4  shows  curves  of  expenditure  of  chlorine  for  con¬ 
ditions  of  temperature  600°K  and  initial  concentration  ratio 
Cl2  ■*  CgHjj,  =  1  :  50,  where  constant  is  computed  by  the  fon 
mula  of  Arrenius  and  not  by  formula  (to).  The  values  38  and 
41  kilocal.  per  mole  {source  203  were  taken  for  the  energy  of 
activation  of  disassociation.  For  an  energy  of  activation  of 
38  kilocal.  per  mole,  81$  of  the  original  chlorine  was  left 
after  10-*sec.;  for  ac  energy  of  activation  of  41  kilocal.  per 
mole,  92.5$  of  the  chlorine  was  left  after  10“+sec. 


5a,l.  Of*  tun  fcmiQ) 


Figure  4.  Curves  of  expen¬ 
diture  of  chlorine  for  E  =  4l(l) 
and  E  =  38  (2)  kilocal. /mole. 
Key: 

1.  [Cl^jin  volumetric  shares 

2.  t  in  seconds. 
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According  to  data  from  source  [21J  ,  for  a  temperature  of 
593°K  and  a  speed  of  supply  of  the  initial  mixture  of  O.323 
moles/liter*sec. ,  the  reaction  time  is  8,10“'isec.  and  84.5$  of 
the  original  amount  of  chlorine  was  left.  This  is  in  satis¬ 
factory  agreement  with  our  calculation.  At  the  same  time  the 
appearance  of  vinylchloride ,  which  was  obtained  as  a  result  of 
the  calcu1  ation,  for  the  same  conditions  already  listed  was 
considerably  less  than  that  experimentally  obtained  by  the 
authors  of  source  .  Evidently  this  is  connected  with  the 

fact  that  the  reaction  process  used  by  us  has  only  one  route 
for  the  formation  of  vinylchloride — reaction  number  13  in 
equation  set  (11) — which  is  insufficient  to  describe  the  given 
process.  Actually,  as  was  already  shown,  even  increasing  the 
magnitude  of  the  constant  of  this  reaction  (kjo)  can  not  in¬ 
sure  an  appearance  of  vinylchloride  comparable Jto  that  ob¬ 
tained  experimentally. 


'T’hls  important  result,  obtained  by  means  of 
computation,  supports  the  well-known  supposition 
"P'n  of  vinylchloride  proceeds  concurrently  by 


a  numerical 
that  the  for- 
two  routes 
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(equations  9  and.  10).  In  other  words  the  formation  of  a  cer¬ 
tain  part  of  the  vinylchloride  is  connected  with  the  reaction 
of  direct  replacement  (equation  10);  a  calculation  for  this 
reaction  is  as  yet  impossible  because  of  the  lack  of  data  for 
the  constants.  Thus  the  question  of  the  relationship  between 
the  reaction  t>f  combination  and  the  reaction  of  replacement  in 
the  process  of  the  chlorination  of  ethylene  for  various  tem¬ 
perature  ranges  needs  further  investigation. 

There  are  also  a  number  of  experimental  results  of  a 
qualitative  nature  from  source  T21J  which  agree  with  the  results 
of  our  calculations: 

1)  Under  the  conditions  of  experiments  described  in 
source  (2i]  HC1  and  C2H-C1  are  formed  in  equi -molecular  amounts. 
The  very  same  conclusion  follows  from  our  calculations.  (For 
example,  see  Figure  la.) 

2)  The  replacement  process  has  a  certain  period  of  induc¬ 
tion  that  decreases  with  an  increase  in  temperature.  A  similar 
dependence  on  temperature  of  the  period  of  induction  was  ob¬ 
tained  as  a  result  of  the  calculations.  (See  Figure  la,  lb). 

3)  With  an  increase  in  time  the  presence  of  products  of 
combination  reaches  a  maximum  and  then  falls  during  the  contin¬ 
ued  development  of  processes  of  replacement.  Such  behavior  of 
o  concentrations  of  products  of  combination,  and  replacement  was 
found  by  means  of  the  calculation.  (See  Figure  la,  lb.) 

4)  In  the  studied  interval  from  600  to  660°K,  the  appear¬ 
ance  of  products  of  replacement  continuously  increases;  the 
appearance  of  products  of  combination  continuously  decreases. 

Our  calculation  results  also  showed  that  the  appearance  of 
products  of  replacement  increases  with  an  increase  of  tempera¬ 
ture,  while  those  of  combination  decrease.  (See  Figure  la, 
lb.) 

5)  In  contradiction  to  the  conclusions  of  source  [2] ,  an 
increased  chlorine  content  in  the  initial  mixture  does  not 
facilitate  a  furtherance  of  the  replacement  process  but  rather 
only  the  formation  of  HC1  and  products  of  intensive  chlorina¬ 
tion.  Our  calculations  also  lead  to  a  similar  conclusion. 

(See  Figure  2. ) 

Figure  5  shows  curves  obtained  for  temperatures  of  600° 
and  800®K,  given  a  ratio  of  CI2  :  C2H4  =  1  :  10  for  the  con¬ 
stant  of  disassociation  of  chlorine;  this  constant  is  1000 
times  larger  than  the  one  computed  by  formula  (16)  and  cor¬ 
responds  to  an  energy  of  activation  of  39.7  kilocal./mole  at 
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600°K  and  3?  kilocal./mole  at  800°K. 


Figure  5.  The  dependence 
on  time  of  concentrations  of 
chlorine,  ethylene,  and  diw 
chlor-ethane,  under  condi¬ 
tions  where  k  =  1000 


T  =  600qK{ continuous  lines) 
T  =  800  K  (dashed  lines) 

1)  £12 »  3L10-Ifc-Ha  , 

5)  io)  jcfHjcrj 

1)  C.  in  volumetric  shares 

2)  tj  in  seconds 


Calculations  were  carried  out  in  which  only  k  _  was 
changed;  all  other  ki  were  held  constant.  Or,  all**^  were 
permitted  to  change  and  k^t;  was  held  constant. 

These  computations  were  carried  out  for  the  very  same 
initial  ratio  Cl2  :  C-H/,  =  1  :  10  for  the  following  cases: 

1)  all  constants  (including  k 1c  calculated  for  600°K;  2)  k15 
calculated  for  600°K,  all  other  constants  calculated  for  8u0°K; 

3)  all  constants  of  reaction  rates  calculated  for  800°K.  Here 
it  turned  out  that  the  time  required  for  expenditure  of  chlo¬ 
rine  sharply  decreases  with  an  increase  in  the  constant  of  re¬ 
action  rate  of  the  disassociation  of  chlorine.  Thus,,  for  in¬ 
stance,  for  15#  expenditure  of  chlorine  the  time  required  in 
the  first  case  was  about  10"1  sec.,  while  in  the  third  case, 
about  10“3  sec. 

Thus,  in  the  process  of  the  chlorination  of  ethylene  ac¬ 
cording  to  reaction  process  (11),  the  value  of  the  constant  of 
reaction  rate  for  the  disassociation  of  chlorine  has  the  great¬ 
est  influence  on  the  process.  Therefore  even  a  relatively 
small  variation  in  the  values  of  the  energy  of  activation  or 
in  the  factor  before  the  exponential  of  the  constant  of  reac¬ 
tion  rate  of  the  disassociation  of  chlorine  may  substantially 
distort  the  results  of  kinetics  calculations,  while  at  the 
same  time,  such  a  variation  in  the  values  of  the  other  con¬ 
stants  of  reaction  rates  evidently  does  not  have  such  an  influ¬ 
ence.  '  Because  of  this  as  accurate  an  experimental  determina¬ 
tion  of  the  constant  k.  as  possible  is  desirable,  as  well  as 
a  study  of  the  conditions  of  the  experiment  on  its  value  (the 
ratio  of  surface  to  volume,  the  nature  of  the  walls  of  the 
reaction  chamber,  etc.*) 
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Til .  A  calculation  of  the  kinetics  of  the  oxidation  of  methane 
and  a  comparison  of  the  results  with  experimental  data. 


The  reaction  of  the  oxidation  of  methane  may  be  consi¬ 
dered  as  a  model  of  a  large  class  of  complex  linked  processes 
of  the  oxidation  of  hydrocarbons .  (Source  [223). 


The  most  valuable  products  of  the  oxidation  of  methane  are 
formaldehyde,  if  the  oxidation  is  conducted  at  pressures  not 
exceeding  atmospheric  pressure,  and  methanol,  if  the  reaction 
occurs  at  high  pressures  (about  100  atmospheres).  Formaldehyde 
and  methanol  are  intermediate  products  of  the  oxidation  reac¬ 
tion;  the  final  products  of  this  reaction  are,  of  course,  car¬ 
bonic  acid  and  water. 


In  many  works  devoted  to  the  study  of  the  process  of  the 
oxidation  of  methane,  relatively  small  production  of  formalde¬ 
hyde  has  been  achieved  for  a  single  passage  of  the  raw  material 
through  the  reactor.  (Sources  (23-28.3).  It  is  well  known  that 
when  this  reaction  is  carried  out  at  relatively  low  temperatures 
and  for  a  short  period  of  contact,  the  methane  oxidizes  practi¬ 
cally  totally  into  formaldehyde.  Flowever  under  these  condi¬ 
tions  the  quantity  of  oxidized  methane  is  small  and  consequently 
the  absolute  production  of  formaldehyde  is  small.  (Source  £24}) 
Raising  the  temperature  or  increasing  the  contact  time,  leads 
to  an  increase  in  the  quantity  of  oxidezed  methane,  but  not  to 
the  increase  in  the  production  of  formaldehyde,  which  itself 
under  these  conditions  oxidizes  easily. 


At  the  present  time  the  best  results  in  obtaining  formal¬ 
dehyde  by  the  direct  oxidation  of  methane  (natural  gas)  by  oxy¬ 
gen  in  the  air  are  achieved  by  the  method  advance  in  source  {25J 
The  basis  of  this  method  is  the  catalytic  action  of  potassium 
tetraborate  on  the  oxidation  of  methane,  and  the  braking  action 
on  the  oxidation  of  formaldehyde  caused  by  a  nozzle  introduced 
into  the  reaction  chamber.  Also  nitrous  gases  are  added  to 
the  reacting  mixture;  these  catalyze  the  oxidation  of  methane. 


The  reacting  mixture  is  composed  of  about  33, T  methane  and  about 
6 T*  air.  The  process  is  carried  out  at  a  temperature  of  about 
1000°Y.  The  time  the  mixture  remains  in  the  reaction  chamber, 
'^cont »  about  0.14-0.16  sec.  The  concentration  of  nitric 
oxides  is  about  10. 6i  by  volume. 


The  production  of  formaldehyde  for  one  pass  through  the 


reaction  chamber  is  about  2.4-2. 6%  of  the  methane  by  volume; 
for  circulation  six  times  the  production  is  about  10.6$  by 
volume. 

Attempts  to  analyze  the  kinetics  of  this  interesting  re¬ 
action  have  been  made  in  source  fy]  for  a  solution  of  the  prob¬ 
lem  of  a  one -dimensional  combustion  of  methane.  The  basis  of 
these  attempts  is  the  method,  proposed  in  source  t33 ,  of  the 
linearization  of  a  set  of  equations  which  describes  the  oxida¬ 
tion  of  methane.  It  is  evident  that  without  obtaining  an  evalu¬ 
ation  of  the  accuracy  of  the  solution  when  linearization  of 
the  equation  set  is  performed,  this  method  can  not  be  used. 

(See  Section  I.)  Using  the  method  developed  in  Section  I  we 
have  solved  on  a  computer  a  set  of  differential  equations  de¬ 
scribing  the  process  of  the  oxidation  of  methane. 

1-  Setting  up  the  problem  and  the  methodology  of  the 
computation. 

A  numerical  solution  is  sought  for  the  set  of  differen¬ 
tial  equations  describing  the  kinetics  of  the  gaseous,  iso¬ 
thermal  oxidation  of  methane  by  a  mechanism  that  shall  be 
examined  below..  In  Section  I  can  be  found  the  detailed  expo¬ 
sition  of  the  algorithm  of  numerical  integration,  to  a  desired 
degree  of  accuracy,  developed  for  complex  chemical  reactions 
for  any  proportions  of  constants  of  reaction  rates  in  the  equa¬ 
tion  set  of  the  kinetics,  the  algorithm  which  we  are  using  here. 

The  chemical  model  used  in  this  article  for  the  process  of 
the  oxidation  of  methane  is  based  on  the  radical-chain  model 
of  the  process  proposed  by  N.  N.  Semenovyy.  (See  Table  2 .*) 


•K-T'he  authors  discussed  the  mechanism  used  here  in  detail  with 
A.  3.  Malvandyan,  V,  Azatyan,  V.  I.  Vedeneyev,  and  S.  S.  Pol¬ 
yak,  and  express  their  deep  appreciation  to  them. 
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(p* 

n.n 

•  CS^eaNtnut 

^kencTHTa  evopoem  * 

1 

CII,  +  0^  C!I,  +  HOj  * 

10“  exp  (-55000//?T) 

2 

CH,+  oicH,00 

10» 

S 

Clt*00  —  C1I,0  +  OH 

10«*exp(-20000/i?r) 

4 

CH4  +  OH-  CH,  +  H,0 

10"  exp(— 8500/j?r) 

S 

CH,0  +  0II-H,0+HC0 

10“  exp  (— 6300//IT) 

< 

CH,0  +  0,-HCO  +  IIOi 

10“  exp  (— 32000/.R7') 

7 

HCO  +  Or*  CO  +  HO, 

10* 

8 

CH4+H0*-*H,0,  +  CH, 

10“  exp  (—19750 /ItT) 

9 

CH,0  +  HO,-  H,0,  +  HCO 

i0“*xp(-8500/*T) 

19 

CO  -f-  OH—  CO,  +  H 

k,i 

10“  exo  (— 7000/717*) 

11 

ch4+h-ch,  +  h. 

1,5.10“  exp(-11200//?D 

12 

CH,0  4-  H-  HCO  -f  H, 

4-10*  exp  (-2000/  RT) 

13 

CH,  +  CHr-*“c,H, 

10“  .  . 

14 

HOj  4*  HOf— ♦  HjOi  -J"  Ot 

lu 

10“ 

15 

OH  +  HOr*  H*0  +  Oi 

k«t 

10“ 

16 

H  +  HO,-  H,  +  0, 

10“ 

17 

ch,  +  ho£.’ch4  +  o, 

io» 

3)  *  Bee  KOKCTSitTW.  M  notnioieHnex  k„ 

bxciot  pasMepaoeni  a'hm.-c**; 

kt  —  «neer  pajuepHOCTb  eer**.  - 

• 

Table  2.  The  mechanism  of  the  oxidation  of  methane* 

Key: 

1.  Iluniber 

2.  Reaction 

2.  Constant  of  rate  of  reaction  rate. 

h,  All  constants  with  the  exception  of  have  the 
dimensions  liler/mole* sec.  The  constant  has 
the  dimension  sec.-’. 

As  can  be  seen,  reaction  .71  is  a  reaction  representing 
the  start  of  the  chain;  react icnc  72-5  and  7-12  are  a  continu¬ 
ation  of  the  chain;  reaction  ,76  represents  the  branching  of 
the  chai* ;  reactions  #13-17  represent  the  mass  breakage  of 
chains  (in  which  reactions  .713  and  represent  a  quadratic 
breakage).  Tn  the  same  Table  2  are  shown  the  constants  of 
reaction  rates  corresponding  to  the  various  reactions.  Most 
.~>e  the  values  for  the  constants  lasted  were  determined  in  a 


temperature  range  not  exceeding  about  1000°K.  The  usual  values 
were  used  for  all  recombinational  chain  breakages. 


Using  the  method  of  quasi-stationary  concentrations,  from 
the  set  of  reaction  processes  and  generally  coinciding  with  it, 
it  is  possible  to  extract  expressions  for  the  change  over  time 
of  the  reaction  (source  X2-* 0).  This  permits  us  to  establish 
that  the  maximum  concentration  of  formaldehyde  depends  only  on 
the  temperature  and  is  determined  by  the  following  relationship: 

fHCHOUj  -  exp  (-8S00/HT)  fCHJ,  (33) 


This  is  close  to  the  result  obtained  experimentally  in  source 

[291: 

u  ••  -  ■  (34) 

'  IHCHOU*  exp  (-10000/nr)  [CH41. 


It  is  essential  to  note  that  this  relationship  is  ob¬ 
tained  for  the  case  where  the  intermediate  product  is  formed 
and  disforms  itself  by  means  of  a  chain  reaction.  It  should 
also  be  pointed  out  that  the  experimental  results  supporting 
the  validity  of  relationship  (3*0  were  obtained  for  tempera¬ 
tures  not  higher  than  900°K  (source  QQ). 

?or  the  computations  we  examined  15  components,  for 
whose  concentrations  we  have  adopted  the  following  designa¬ 
tions: 

Coeflimenne  CP  CH4  0,  CH,  HO,  CH.OO  CH,0  OH 
KonnenTpanmr^  c,  c,  c,  q  «,  c,  «, 

CoeflnBDHno®  H:0  HCO  CO  H,0,  H,  CO,  H  C,H, 

KoBuciiTpaun^  r»  c,  q#  c„ "  c,4  q,  q,  q. 

Key: 

1.  radical  3.  radical 

2.  concentration  h.  concentration 


Corresponding  to  the  chemical  model  chosen,  a  set  of  15 
non-linear,  ordinary  differential  equations  (35/  was  set  up: 
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iji  as  « —  frjCjC, —  kgCtCt jCiCt 

~~  =  —  *i<y»  —  A-,C,C3  —  Vi'i  —  AfjCjC,  +  fcl4cj  4- 

4*  4"  Aa«c4c44  4-  kiiC^Ct 

=  —  AjCjCj  —  2Auc’  —  ktfCjCi  +  Ar4c4c,  4-  £4fiiC»  + 

4*  A*Cj<r4  -f-  Afncicu 

■y-  —  - — ■  ktfiyCi  — -  A,c4c4  — ■  2Ar44c*  A44c4cj  —  A44c4£44  — — 

- —  AjfCjC4  +  A;CaCj  -)-  AjCjCa  -J-  k^cfy 

~-  =s  —  AjC4  +  AjCjCj 

Jgjp  s=  — ■  kffitfii  — *  AfCjCg — A4£4c4  A|jC(Cj4  4  Ms 

-^-  =  — A4cjc7  —  kifiiCy  —  A10CjCio  —  Acj4c4Cj  -}-  Aj£j 

Hr  —  At4CiCj  +  AsC|C7  -f-  knCtc7 

■4r  =  —  AjCjc,  +  A*c*cj  -f  AjfjCj  4-  A»c4e«  +  Arji£4cl4 
o» 

-jr  =  —  AjoCjCxo  4“  A'jCjCi 

^  =  W4  +  ^«  +  M 

=  Allc1cu  4-  A1]c«cu  4-  A14c4cl4 

*=  Ajo«jC10 

=  —  A11e1Ci4  —  kltctCu  —  4*  AqoC?cio 

ol 
*u 

1C 


+ 


(35) 


Xn  the  case  where  the  reaction  proceeds  in  a  current 
flow,  the  proper  equations  of  hydrodynamics  should  be  added 
to  a  differential  equation  set  (35)  and  the  aggregate  set 
should  be  solved  together,  as  was  done,  for  instance,  in 
source  £l5j.  df  each  of  the  equations  of  the  set  is  pre¬ 
sented  in  the  form 


jfy _ 


dt 


«'(<.  y)y  +  b'  V), 


then,  as  was  shown  in  Section  1,3,  when  the  following  con¬ 
dition  holds: 


a»> 


|!ne| 


whore  "T  is  the  reaction  tine. and  L  the  required  accuracy, 
several  components  will  enter  the  so-called  quasi-stationary 
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state  within  a  time  short  by  comparison  with  the  reaction  time. 
For  this  reaction  which  we  are  examining— the  oxidation  of 
methane — such  components  are  c*  (since,  for  instance,  when  the 
temperature  is  1000°K  ai  -  0.454* 10y)  and  also  under  certain 
circumstances  and  Cgtin  particular,  when  c^  and  c2  are 
sufficiently  large  at  the  start  of  the  reaction).  Here  Ct 
enters  the  quasi-stationary  state  at  a  time  of  the  order  of 
10“6 .  while  c 7  and  Cq— if  the  conditions  concerning  c^(0)  and 
02(0)  are  fulf illed— enter  the  quasi-stationary  state  at  a  time 
one  or  two  orders  of  magnitude  later.  As  was  shown  earlier  in 
Section  1,3,  a  solution  of  this  set  of  equations  using  usual 
methods  is  practically  impossible,  since  the  indicated  equa¬ 
tions,  for  entry  of  the  corresponding  concentrations  into  the 
quasi-stationary  state,  yield  oscillating,  unstable  solutions 
in  the  numerical  integration  process.  The  algorithm  used  by 
us  enables  us  to  remove  this  difficulty  entirely. 


It  is  important  that  in  this  case  the  equations  for  the 
determination  of  c c,  Co,  and  Cq  are  independent,  of  each  other; 
this  considerably  simplifies  the  computation.  The  difference 
between  the  value  of  concentration  Cr  from  its  quasi-stationary 
value  becomes  negligibly  small  when  \lnt\  /k„.  Moreover, 
inasmuch  as  and  c  take  part  in  second-order  reactions  (and 
consequently ,  at  may^decrease  as  a  function  of  the  rates  of 
disassociation  6.,  Cg,  c.0>  c^,;  likewise  at  as  a  function  of 
the  rate  of  disassociation  of^Cp),  they  may  with  the  passage 
of  time  leave  the  quasi-stationary  state.  Hence,  in  order  to 
insure  the  desired  degree  of  accuracy  in  the  computation  pro¬ 
cess,  it  is  essential  to  check  for  the  satisfying  of  the  con¬ 


ditions  for  the  quasi-stationary  condition  for  c  and  c 


9' 


The  numerical  integration  of  this  set  of  equations  was 
done  by  the  method  of  Rungs— Kutt  with  automatic  selection  of 
the  interval  of  integration  and  a  relative  error  of  from  10“^ 
to010“-6  Initial  concentrations  of  methane  and  oxygen 
(C^  :  C2  *  O.29  :  0.71)  are  taken  to  be  the  initial  conditions 
of  the  set  of  differential  equations  (35)  that  is  under  exami- 
naf-i  on. 


2.  discussion  of  the  results  of  the  numerical  integra¬ 
tion  of  set  (35)  and  a  comparison  of  them  with  experimental 
data. 

2a)  The  functions  c.  =  c.(T),  c.  =  c.(t).  dc./dt  =  f(t). 

J  J  J  J  1 

The  numerical  integration  of  equation  set  (35)  was  car- 
r i, 0 A  out  for  the  temperature  range  300-1700°K,  As  a  result, 
the  dependence  on  time  of  the  concentrations  of  all  substances 
taking  part  in  the  reaction  was  obtained  for  various  tempera- 
in  this  range,  as  well  as  curves  of  their  expenditure 
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and  formation.  The  dependence  on  temperature  of^all  concen¬ 
trations  that  was  thus  found,  for  the  instant  10"*4,  sec.  is 
shown  in  Figures  6 a  and  6b.  Moreover,  in  order  to  clarify 
the  influence  of  chain  breakages  (see  equations  # 13-17  in  Table 
2)  on  the  process,  calculations  were  carried  out  under  the  same 
conditions  for  the  reaction  process  1-12. 


tUr  ICH‘,: '  “  r°,1:  *  -,0‘ [CH,03: W— totcoj;  is — 10*  [H,l;  n-to’ 
ICOtl* 

wiToMH?  14-To!»°iShJ  ’  ,~10'  ICH*00,:  ,-,8‘ l0Bl*  *-10'  {HC°i 


Figures  6a,  6b 

Tine  dependence  on  temperature  of  concentrations  of 
substances  taking  part  in  the  oxidation  of  methane. 

Kcy:r 

1.  (jpHhj  in  volumetric  shares. 

2.  c..( j  =  2)  in  volumetric  shares. 

3.  cj  in  volumetric  shares. 


Figures  ?a,  7b  and  8a  and  8b  show  in  semi-logarithmic  scale 
the  computation  results  at  1300°K  for  the  abbreviated  (not 
taking  into  account  chain  breakages)  and  full  reaction  process, 
and  Figures  9a  and  9b  show  the  same  curves  at  300°K  for  the 
full  reaction  process. 

It  follows  from  an  examination  of  Figures  7a  and  7b  that 
the  curves  for  concentrations  of  Cp,  HOo,  H,  iUC^HgO,  H2  and 
CCo  reach  saturation  at  a  time  of  the  order  of^lO-3  sec. 
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a)  1  —  [CH.l;  F-CO.l;  J— lO'fCH,];:  <— 10*  tHO.J:  «— 10*  [CH.Ol;  7-10*[OH \jl— 

(H1O1I. 

'  i-io*  tCH.OOl;  t  —  IH.OI;  »-i0<  iHCO];  It  -  fCO];  It  —  [H,];  13  -  ICO,); 

H-10*  Ill] 

Figures  7a,  7b.  The  dependence  on  time  of  the  concentra¬ 
tions  for  the  reaction  set  from  Table  1  (without  reactions 

of  chain  breakage) 

Key: 

1.  c .  in  volumetric  shares. 


)  <-  ID*  fHOj;  f— 10*  [CH,0];  7-10*  [OH];  F-10*  (HCOJ;  It  -  [H,l: 

M  —I  UCOi);  14—10'  [Hi 


igures  3a,  8b 


The  dependence  on  time  of  the  concen- 


t rat ions  for  the  complete  set  of  reactions  # 1-17  from  Table 
2  (including  reactions  of  chain  breakage). 

Key : 

1.  c.  in  volumetric  shares. 

2.  in  volumetric  shares. 

3.  tJin  seconds. 

4.  t  in  seconds. 

Tn  the  case  where  the  complete  set  of  reactions  is  examined 
(see-  Figures  8a,  8b),  as  could  have  been  expected,  the  curves 
for  the  concentrations  of  the  radicals  (such  as  HOp  and  H) 
pass  through  a  maximum,  while  the  concentrations  or  f^O,  HLO2 
[ip  and  others  (which  are  the  products  of  the  reactions  of  * 
chain  breakage  continue  to  increase.  For  the  same  tempera¬ 
tures  and  for  the  set  of  reactions  // 1-12  (see  Table  2),  cal¬ 
culations  were  carried  out  for  a  decrease  in  constant  by 
a  factor  of  1000,  which  corresponds  to  an  increase  in  the 
energy  of  activation  of,  for  example,  14  kilocal.  at  1000°K 
and  almost  18  kilocal.  at  1300°K.  Here  it  would  seem  that 
such  a  decrease  in  the  value  of  k„  has  practically  no  effect 
at  all  on  the  process.  This  situation  agrees  with  N.N. 
Semenov's  proposition  that  the  process  of  the  formation  of 
formaldehyde  is  determined  by  the  reaction 

CH.  + 0^01,0  + OH, 

since  reaction  #3  follows  quite  quickly  after  reaction  If 2 
(See  Table  2)  (Source  DO)* 


As  could  have  been  expected,  with  an  increase  in, tempera 
turc  the  curves,  while  maintaining  their  likenesses,  are  dis¬ 
placed  ("compressed,"  and  "stretched  out")  along  the  time  axi 
(dee  Figures  8a,  Sb,  9a,  9b.) 

-  Sa 

aflW/w  Do/iufe  .  Cj,  oOivum.  flt/miQ 

*•  $r~ - - 1 


to”  io  -  m  •  io''  - 

.  •  -  t  czk®  f-rvuyff 

»_7o S;,  *  -  t0l,;  4  -10’  lCH']=  '-*»■  tHO.);  *_!0»  IHC01;  iHl; 

tO>  [H.OJ;  Jo  10*  [CO];  Jl~ 


Figures  9a,  9b.  The  dependence  on  time  of  the  concen¬ 
trations  for  the  full  set  of  reactions. 

Key: 

1.  a  .  in  volumetric  shares. 

9.  in  volumetric  shares. 

3.  fcJin  seconds. 

4.  t  in  seconds. 


Figures  10  and  11  show  the  dependence  on  time  of  the  ex¬ 
penditure  of  CH^  and  0„  for  temperatures  of  1000°K  and  1500°K. 
For  purposes  of  comparison  the  curves  of  dependence  on  time  of 
the  concentration  of  formaldehyde  corresponding  to  these  tem¬ 
peratures  has  been  added  to  these  figures.  From  an  examination 
of  Figures  10  and  11  it  can  be  seen  that  with  an  increase  in 
temperature  the  maxima  of  the  expenditure  curves  of  CH^  and  0^ 
get  closer  (in  time)  to  the  maximum  of  the  curve  for  the  con- 


Figure  10. .  The  dependence 
on  time  of  the  rate  of  ex¬ 
penditure  of  1)  methane, 

2)  oxygen,  and  3)  the  con¬ 
centration  of  formaldehyde 
for  a  complete  set  of  reac¬ 
tions  at  1000°K. 

Key: 

1.  t  in  seconds. 


Figure  11.  The  dependence 
on  time  of  the  rate  of  ex¬ 
penditure  of  1)  methane, 

2 )  oxygen,  and  3)  the  con¬ 
centration  of  formaldehyde 
*  for  a  complete  set  of  reac¬ 
tions  at  1500°K. 

Key: 

1.  t  in  seconds. 
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Figure  12  shows  the  dependence  on  time  of  the  rate  of  for¬ 
mation  of  CH„0  for  different  temperatures.  With  an  increase 
in  temperature  the  curves  are  displaced  leftward  along  the 
time  axis  in  the  direction  of  smaller  times;  the  maximum 
values  of  the  curves  increase;  but  their  forms  do  not  change. 


*  it 


Figure  12 .  Dependence 
on  time  of  the  rate  of 
formation  of  formalde¬ 
hyde  for  different  teim 
peratures. 

Key: 

1.  t  in  seconds. 


/  —  T  »  isoonc:  a  «  !0-»;  *  —  r  ->  i300*k;a  -  to—*:  *  —  T  *= 

~  1000’K.  A  =■  10-';  4  —  T  —  SOO-K.  A  —  1« 

In  order  to  study  the  change  in  the  mechanism  of  the 
oxidation  of  methane  in  the  temperature  range  under  considera¬ 
tion,  the  length  of  the  chains  was  calculated.  The  numerical 
values  for  the  calculation  were  taken  at  points  where  the  con¬ 
centration  of  formaldehyde  was  90 #  of  its  maximum. 

The  results  of  the  computation  were  as  follows: 

T,  *R  .......  1000  1300  1500  1700 

(p  v  (nnnna  rteim) .  .  .  120  77  17  0,41 

Key : 

1.  v(the  length  of  the  chain) 

It  is  evident  from  these  data  that  the  chain  mechanism 
for  the  oxidation  of  methane  occurs  only  up  to  temperatures  of 
about  1600~1700°K,  in  which  range  equation  (33),  obtained  on 
the  basis  of  a  chain  process  of  reactions,  is  already  no  longer 
in  effect.  We  have  no  data  of  any  kind  at  our  disposal  con¬ 
cerning  the  type  of  mechanism  for  the  oxidation  of  methane  at 
higher  temperatures.  It  may  be  supposed  with  some  degree  of 
certainty  that  it  will  be  of  the  free-radical  type. 

2b)  The  energy  of  activation  of  the  oxidation  of  methane. 

As  is  well  known,  the  energy  of  activation  of  the  oxida¬ 
tion  of  methane  decreases  with  an  increase  in  temperature.  The 


corresponding  data  are  presented  in  Table  3.  in  which  the  con¬ 
ditions  of  the  experiment  are  also  indicated,  since  the  numeri¬ 
cal  values  obtained  also  depend  on  the  condition  of  the  surface 
of  the  reaction  chamber  and  the  nature  of  the  mixture. 
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CH4  +  20. 

Coeya  lift 

o8pa8otau  HP 

32(773)— 25(798) 

TaS.i.  4CJ 

748-798 

6,63 

JCH4+2O, 

Cocyff  jTJv 

•crapuft*  ViV 

9S  (748)— 65  (798) 

898—948 

873—948 

0,OG6 

0,13 

CH«  +  O, 

2CH4h30. 

Coeya  oSpafio- 
in  HP 

Coeya  nporper 

43 

78 

122, 

tb6.i.  47] 

848—935,5 

0.4 

2CII4+30, 

Coeya  nqijpuT 
pbo 

57—43 

1 57 

873-948 

0,2 

2CII4  -f-  Oi 

Coeya  nporper 

873-935,5 

0,3 

2CH4+30. 

Coeya  noKofcrr 
p»o  qg 

70—55 

1050- 

70%  Ar  -f- 

21,5-33,8 

J32J 

2100 

9%CH4  + 
21%  O, 

‘ 

1420- 

4,  <P< 

70%  Ar  + 

49,5 

[32] 

2000 

>71 

18%CH4  + 
12%  O, 

1100- 

2000 

Eojiee 

Ncrftiia  b  mb- 

TB0O-UH&1O* 
POflHOW  CMCCM 

52,9 

[33] 

700—1700 

Mexee  50% 

MfrTdun  ■  ue- 
TBHQ-KHCyiO' 
POflHOft  CMCCM 

*(3) 

20,6 

[33] 

•  B  CKoGKa^jrKaaanw  teKnepaTypw,  jouirt  kotopux  fiO-iyieim  oimenrtn  aneprnH 

AKTHBaUHK. 


Tabic  3*  '"he  dependence  of  the  energy  of  activation  upon 
the  conditions  of  the  experiment. 


Key: 

i . 

9  > 

3! 

/i. 

n; 

r.\ 

i  • 

o 

o 

to! 


temperature  range,  °K. 
conditions  of  the  experiment. ' 
t>  in  atmospheres, 
nature  of  the  mixture, 
notes . 

E  in  kilocal./mole. 

Source  of  data. 

Reaction  ohamber  treated  vrith  hydrofluoric  acid. 
Reaction  chamber  "old." 

Reaction  chamber  treated. with  hydrofluoric  acid. 
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11.  Reaction  chamber  "old." 

12.  Reaction  chamber  treated  with  hydrofluoric  acid. 

13.  Reaction  chamber  heated  up, 

14.  Reaction  chamber  covered  with  lead  oxide. 

15.  Reaction  chamber  heated  up. 

16.  Reaction  chamber  covered  with  lead  oxide. 

17.  table  46,  table  47,  etc. 

18.  More  than  54#  methane  in  a  methane- oxygen  mixture. 

19.  Less  than  50#  methane  in  a  me thane -oxygen  mixture. 

20.  *  The  numbers  in  parentheses  are  the  temperatures  for 
which  the  values  of  the  energy  of  activation  were 
obtained.  . 


For  all  of  the  picture's  complexity,  it  can  be  seen  that  the 
energy  of  activation  decreases  with  an  increase  in  temperature 
(and  also,  evidently,  with  an  increase  in  pressure). 


Source  C22l  offers  several  reflections  on  this  subject. 
At  any  rate,  the  indicated  tendency  is  distinctly  noticeable 
in  all  of  the  experiments. 


In  this  connection  let  us  examine  the  energies  of  activa¬ 
tion  of  the  oxidation  of  methane  shown  in  Table  4,  which  were 
derived  from  the  relationships : 


In 

In 

In 

In 


<f  (CH« 


dt 

rfro.11 

rf*  ItTllX 

rf[C0]| 


dt  jrt 

rf[H.0| 


dt 


r'(i) 

-'(« 

-'W 

Mr) 


V). 

(*>. 

(5). 

<*). 


Burepun 


E,  kkqa/moa b 


© 


tenrirjaTyp. 

I 

J 

4 

1200-1300 

1300-1400 

1500—1600 

25,0 

25,0 

18,6 

21,2 

23,62 

18,6 

30,0 

26,6 

18,2 

32,8 

28,0 

17,0 

Key: 

1.  Temperature  range, 

2.  E.  in  kilocal./mole. 


K. 


Table  4,  Various  values  for  the  energy  of  activation  of 
the  oxidation  of  methane. 


It  can  be  seen  front  Table  4  that  the  tendency  (noted  in 
the  experiments  listed  in  Table  3)  for  the  energy  of  activation 
of  the  oxidation  of  methane  to  decrease  for  an  increase  in  tem¬ 
perature — that  this  tendency  derives  from  the  mechanism  shown 
in  Table  2,  because  the  data  of  Table  4  are  the  mathematical 
consequence  of  a  computation  of  that  mechanism. 

The  values  for  the  energy  of  activation  calculated  by 
different  methods  agree  well  with  each  other,  with  the  excep¬ 
tion  of  the  values  determined  in  the  1200-1300°K  range  using 
the  expenditure  of  CH^  and  0?,  on  the  one  hand,  and  using  the 
formation  of  CO  and  HgO,  on  the  other  hand.  The  reason  for 
this  discrepancy  is  not  clear .although  the  problem  was  noted 
in  several  sources  (source  (33 )  • 

2c  )  Comparison  of  computation  results  with  experimental 
data. 


An  experimental  investigation  of  the  oxidation  of  methane 
was  conducted  in  a  plasma  stream  of  argon  (sources  C3df  £3if| ) . 
We  shall  compare  these  experimental  data  with  the  results  of 
the  computation.  First  of  all  it  must  be  noted  that  the  results 
of  experiemnts  in  a  stationary  plasma  stream  at  atmospheric 
pressure  for  various  chemical  reactions  agree  well  with  data 
obtained  under  static  conditions  (source  ObJ) ♦  The  error  of 
temperature  measurement  in  sources  £34113 5J  was  i  10  -  12$, 
which  does  not  exceed  the  usual  errors  in  that  temperature 
range;  the  error  in  measuring  the  concentrations  of  most  of 
the  gases  was  +  5,* ,  and  that  in  the  determination  of  CH_0  was 
+  20^-  ■  2 


The 


following  can  be 


seen  from  an  examination  of  Figure 


13: 


1  -  ^jjjjjtC-sMnepKsieiiffii&iiMe  Town);  l-  (CJ-  wfcnepi^rwtffiie  town); 

*  —  '  ^'(X—SKcncpH^^TaJiMme  njMMH) 


Figure  13 
A  comparison 
of  calculated 
curves  with  ex¬ 
perimental  data. 
Key: 

1 .  experimental 
points . 

2.  experimental 
points. 

3.  experimental 
points. 


-  46  - 


1)  the  experimental  points  for  the  temperature  interval 
of  calculation  from  800  to  1700  K  satisfactorily  agree  with 
the  calculated  curves. 

2)  Both  the  calculated  curve  and  the  experimental  data 
for  the  formation  of  formaldehyde  show  a  maximum  approximately 
in  the  1500-1600°K  range.  Thus,  this  important  result  of  the 
calculation  is. fully  confirmed  by  experiment.  It  is  essential 
to  note  that  this  maximum  determines  a  formerly  unknown  fron¬ 
tier  in  the  usefulness  of  equation  (33)»  It  remains  an  open 
question  whether  this  frontier  is  connected  with  a  change  in 
the  mechanism  of  reaction  or  occurs  for  some  other  reason. 

3)  Calling-  attention  to  itself  is  the  fact  that  production 
of  CH?0  (curve  2  in  Figure  13)  grows  considerably  faster  up  to 
the  maximum  point  than  it  decreases  after  the  maximum.  This 

is  easily  seen  from  an  examination  of  the  experimental  data  and 
is  noticeable  on  the  calculated  curve  (although  the  part  of  the 
calculated  curve  after  the  maximum  is  not  large). 

Let  us  now  compare  the  calculated  data  shown  in  Table  4 
with  the  result  of  a  calculation  of  the  energy  of  activation 
found  in  source  D  53  for  pressure  p  =  1  atmosphere  (See  Table 
5). 


£.  HttttAjMOAb 

Lx 

lIiiTopna.i 

TvKiena- 

t 

|  f 

1 

TJFp,  ‘K 

W'liicS? 

mJP 

© 

•HOIItptIM. 

w 

BU<0(0.1. 

CD 

BKcnepHM.'-^ 

1200—1300 

1500—1600 

2T.-10* 

18.6* 10* 

(29,6±«)  10* 
(I0.0±3.3)- »• 

24,2*10’ 
13, C* 10*  j 

_t3,9)IO* 

MP.6±2,t}IO’ 

30.0- 101 

18,2- 10«  j 

<28,5.fc4,S)l0» 

(I8.0±3,l 

Table  Values  of  the  energy  of  activation  determined 
by  1)  the  expenditure  of  CHg ,  2) the  expenditure  of  C0  and 
3)  the  formation  of  CO. 


Fey: 


1. 

t  '-mperature  range, 

5. 

calculated 

2. 

11  in  kilocal  ./mole 

6. 

experimental 

3. 

calculated 

7. 

calculated 

it. 

•  • 

experimental 

T. 

experimental 

If  the  inaccuracies  of  the  experiment  and  of  the  initial 
data  for  the  calculation  are  taken  into  account,  then  the  agree¬ 
ment  between  calculated  and  experimental  values  of  the  energy 
of  activation  must  bo  conceded  to  be  a  good  one.  Thus,  the 
r* v per  1  mental  data  and  the  calculation  according  to  the  mechanism 


^5'^'ps^ivw»*r  "* 


of  Table  2  agree  and  both  indicate  what  has  already  been  men¬ 
tioned  by  many  authors--that  the  energy  of  activation  of  the 
oxidation  of  methane  decreases  with  an  increase  in  temperature. 

In  conclusion  we  shall  make  one  comment  in  the  spirit  of 
the  proposition  developed  later  in  this  book  Cin  a  chapter 
much  later  than  the  chapter  here  being  translated. .TransH. 
about  using  the  computer  to  search  for  the  most  probable  mech¬ 
anism  of  this  or  that  chemical  reaction.  For  the  calculation 
of  the  mechanism  of  reaction  of  the  oxidation  of  methane,  lis¬ 
ted  in  Table  2,  we  have:  a  good  approximation  by  experimental 
points;  the  absence  of  the  influence  of  K_  on  the  oxidation 
process;  the  fulfilling  within  discrete  limits  of  equation  (3*0; 
the  agreement  of  the  calculated  values  for  the  energy  of  acti¬ 
vation  with  experimental  observations ,  and  the  fulfillment  of 
the  experimentally-found  tendency  for  the  energy  of  activation 
to  decrease  with  an  increase  in  temperature.  Although  we  did 
not  examine  other  possible  mechanisms  for  the  oxidation  of 
methane  (this  was  not  the  purpose  of  the  present  article),  a 
review  of  the  results,  which  were  not  earlier  inserted  into 
the  computer  program,  shows  that  the  mechanism  of  Table  2  is, 
at  the  least,  one  of  those  which  correctly  reflect  both  the 
totality  and  the  sequence  of  stages  of  the  reactions  of  the 
real  process  of  oxidation  of  methane. 

IV.  The  influence  of  initial  conditions  on  the  solution 
of  a  set  of  equations  of  chemical  kinetics  and  hydrodynamics, 
(for  example  the  pyrolysis  of  methane  in  a  current  stream). 

An  investigation  of  chemical  processes  (especially  the 
nen- isothermal  ones  in  a  current,  which  are  one  of  the  more 
■widely  encountered  types  in  chemical  technology) ,  demands  a . 
clarif ication  of  the  effect  of  a  change  in  the  initial  condi¬ 
tions  (temperature,  the  ratios  of  concentrations,  etc.)  on 
t>0  solution  of  the  differential  equation  set  that  is  the  model 
of  the  process.  In  this  section,  using  the  pyrolysis  of  methane 
in  a  plasma  stream  as  an  example,  such  an  investigation  is  car¬ 
ried  out  by  means  of  a  numerical  solution  on  an  electronic  digi¬ 
tal  computer.*  In  source  [153,  on  the  mathematical  model  of 
a  plasma-chemical  process  of  the  Conversion  of  methane  to 

*In  an  appendix  Cto  the  book  in  which  this  article  is  found, 
on  p.  233  the  question  of  using  electronic  analog  computers 
for  problems  of  this  type  is  discussed. 


acetylene,  the  following  items  were  taken:  the  initial  tempera¬ 
ture"  T ( 0 ) ,  the  initial  speed  v(0)  of  the  plasma  stream,  the 
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initial  concentration  ci  (0)  of  methane;  and  the  influence  of 
these  items  on  the  maximum  concentration  of  acetylene  o»,(z^), 
the  length  of  the  reaction  chamber  L,  and  other  magnitudes 
was  studied.  As  a  result  it  was  noted  that  the  dependence  of 
Ci(zw),  and  other  variables  on  the  initial  conditions  enables 
us  to  take  a  look  at  the  possibility  of  optimally  controlling 
the  plasma-chemical  process  by  means  of  corresponding  changes 
in  these  [initial]  conditions. 

VJe  regard  the  plasma-chemical  reaction  chamber  as  a  device 
of  automatic  control.  It  is  well  known  that  in  order  to  describe 
a  device  of  automatic  control,  it  is  necessary  to  specify  the 
relationships  between  its  outputs  and  inputs.  These  relation¬ 
ships  may  be  found  for  a  plasma-chemical  reaction  chamber  by 
using  the  mathematical  model  of  source  ^53of  the  process  of 
the  conversion  of  methane  to  acetylene.  In  this  source  was 
examined  the  following  set  of  ordinary,  non-linear  differential 
first-order  equations,  which  to  certain  degrees  of  approxima¬ 
tion  simulate  the  plasma-chemical  reaction  of  the  conversion 
of  methane  to  acetylene  in  a  hydrogen  plasma  stream: 


<r> ~  *■  (*)] 
nr1  *  W)  k' {T)  w  • “  *»  W  *  <*>] 

+J~Mr>c3(2)] 

f  AT  ~  2v  w]  S  ci^  -Jf  +  9  (J)  *  2  Cj  UWj 


*1  _  1 
di  »(z) 


where 


- 


(N  and  M_ are  constants  determined  by  the  initial  condi¬ 
tions.) 


The  numerical  integration  on  a  computer  of  the  equation 
set  (36),  for  a  given  aggregate  of  constants  of  reaction  rates 
k, (T) ,  ko.(T).  and  k$(T) ,  enables  us,  by  changing  the  initial 
conditions  T(0),  v(0),  and  c  (0),  to  establish  a  dependence 
on  T(o),  v(0)f  c,  (0)  of  the  solutions  of  this  equation  set  for 


any  .z.  In  connection  with  this  it  is  natural  to  take  the  mag¬ 
nitudes  T(0),  v(0),  and  'c.  (0)  as  inputs  to  the  plasma-chemical 
reaction  chamber.* 


*We  shall  not  further  be  interested  in  the  dependences  of  the 
outputs  on  input  v(0),  since  it  is  shown  in  source  [l^fjthat  a 
change  in  v(0)  does  not  influence  the  kinetics  of  the  chemical 
reactions,  but  only  changes  the  special  extent  of  the  process. 


The  choice  of  magnitudes  [yariablesj  to  be  taken  as  outputs 
depends  on  the  end- purpose  of  the  process  and  is,  to  a  certain 
extent,  given.  As  quantities  to  be  taken  as  outputs  of  the 
plasma-chemical  reaction  in  this  article  the  following  will  be 
examined,  all  determined  at  the  point  where  the  concentration 
of  acetylene  reaches  a  maximum  (zm): 

a)  the  length  of  the  segment  measured  from  the  end  of  the 
reaction  chamber  z  =  O'  to  z  =  z  .  This  segment,  in  which  the 
reaction  of  the  conversion  of  methane  to  acetylene  takes 
place  for  the  most  part,  is  designated  L  and  is  called  the  ef¬ 
fective  length  of  the  plasma-chemical  reactor,  in  accordance 
with  the  term  used  in  source  ; 

b)  the  maximum  of  the  concentration  of  acetylene,  *0^(2^); 


c)  the  degree  of  decomposition  of  methane. 


* 

Each  of  the  outputs  is  a  function  of  all  of  the  inputs, 
and  also  parametrically  depends  on  the  constants  of  reaction 
rates  ki,  k? ,  and  k~,  which  are  coefficients  in  equations 
set  (367.  That  is. 

£=*iir(0),  cl(0),v(0);  *1.  fc»]  ] 

e"i(2m)=£»  1^(0).  <h(0),  y  (0);  *1,  &s]  [•  (3®) 

*  =**ir(0)A  (0).  y (0);  h,  kt,  k3 ]  J 
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Because  equation  set  (36)  can  not  be  solved  analytically, 
it  is  impossible  to  determine  directly  the  functional  depen¬ 
dences  (3?)  of  the  outputs  using  the  inputs  and  constants  ki, 
kg,  and  ko.  However,  numerical  integration  by  computer  of 
our  equation  set  enables  us  to  solve  the  following  problems: 

1)  The  construction  of  the  approximate  dependence  of  the 
outputs  on  the  inputs  for  a  given  choice  of  the  constants  k. , 
kx  and  k3.  These  dependences  will  be  called  the  static  charac¬ 
teristics  of  the  plasma-chemical  reactor  in  the  discussion  to 
follow. 

2)  The  clarif ication  of  the  degree  to  which  the  static 
characteristics  of  the  plasma-chemical  reactor  change  for  a 
change  in  the  aggregate  of  constants  k>  ,  k-*.,  and  k-a,  and  as  a 
result  of  this,  the  examination  of  the  influence  of  an  inac¬ 
curate  determination  of  the  constants  of  reaction  rates  on 
the  outputs  of  the  reactor.  Differences  in  the  values  of  the 
outputs  for  fixed  inputs,  caused  by  changes  in  one  or  another 
of  the  constants  of  reaction  rates,  may  serve  as  a  measure  of 
the  sensitivity  of  the  mathematical  model  of  the  plasma-chemi¬ 
cal  reactor  to  such  changes  in  these  constants  of  reaction 
rates . 


In  view  of  the  fact  that  under  real  conditions  the  inputs 
undergo  random  perturbations ,  the  problem  arises  of  studying 
the  effect  of  these  perturbations  on  the  outputs  of  the  reactor. 
He  consider  this  problem  as  a  problem  of  the  passage  of  a  ran¬ 
dom  perturbation  through  a  transformer,  the  latter  in  this 
case  being  the  plasma-chemical  reactor.  The  mathematical  expec¬ 
tations  concerning  the  outputs  will  depend  on  the  statistical 
characteristics  of  the  perturbations.  It  is  evident  that  these 
dependences  need  not  be  taken  into  account  when  the  variation 
in  the  value  of  the  outputs  caused  by  an  inexact  determination 
of  the  constants  of  the  reaction  rates  is  considerably  larger 
than  the  changes  in  outputs  caused  by  random  perturbations, 

1,  Setting  up  the  static  characteristics  of  the  plasma- 
chemical  reactor.  ~  ~ 


By  numerically  integrating  the  equation  set  (36)  for  five 
variants  of  each  of  the  inputs  c  (0)  or  T(0)  with  tiie  other 
input  held  constant,  it  is  possible  to  construct  the  graphical 
dependences  of  all  of  the  outputs  L,  c5 ,  z*  and  s_  on  each  of 
the  inputs  c, (0)  and  T(C). 

Figure  l't-a  shows  the  dependence  of  c^(zm),  s_  and  on 
T( 0 )  for  fi#ed  values  of  inputs  <S\  (0)  =  0.95  by  weight  and 
v(0)  -  3*10  cm/sec.  Figure  14b  shows  the  dependence  of 
c  (z  ),  s  and^L  on  c  (0)  for  fixed  values  of  inputs 
T(0)  =  3.5*10J  °K  and  v(0)  =  3-104  cm/sec. 


- 


14  b 


L,C*  Zj  (2f,l  lie.  Buna 


Figures  14a,  l4b.  The  dependence  of  l)  the  maximum  con¬ 
centration  of  acetylene,  2)  the  degree  of  decomposition 

of  methane,  s,  and  3)  the  length  of  the  reactor,  Lj  on  dif¬ 
ferent  ,  inputs'  into  the  plasma-chemical  reactor. 

Key  5^ 

1.  0-3(2*),  by  weight. 

2.  ViM,  by  weight. 

3.  T(0).  10J  °K 

4.  c,(0j,  ty  weight 


Using  the  curves  it  is  possible  to  write  the  approxi¬ 
mate  analytical  dependences  between  the  given  inputs  and  out¬ 
puts,  by  means  of  approximating  these  curves  with  linear  and 
quadratic  functions. 

The  dependences  of  outputs  d(z *) ,  and  s_  on  the  outputs 
T(0)  and  c  (0)  are  shown  by  the  corresponding  formulas  (39-41) 
and  (42-44 ) . 
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(39) 


Mzm)  -  0,125  +  3,5  •  10“4  {T  (0)  -  2000}, 

(2000<7’(0)<3000°K1.* 

c,  (zm)  =  0,476  +  3,44  •10'4{7'(0)  —  3000}  — 1,6-  10-7{T(0) — 3000}*, 

(39a) 

[300p<7’(0)<3500"Kl. 

L  =  67— 9- 10~*  {T  (0)  —  2000}  +  0,44- 10-4  {T  (0)  —  2000}*,  740) 
{2000  <  ?(0)<  3000"  KJ. 

L  =  21—2,6 •  10"*  {T (0)  —  3000}- 0,8 ■  10'* {f(0) — 3000}*,  (40a) 
{3000  <2'(0)<  3500  °KJ. 

*  =  24,7  +  5,83- 10~*  {T  (0)  -  2000},  (41) 

(2000  <r(0)<  3000  °K}. 

a  —  83  +  0,558  ■  10'*  {T  (0)  —  3000}  —0,464  •  10‘4  [T  (0)-  3000}*, 

(41a) 

[3000  <  r(0)<  3500  °K]. 

c*  <*»)  =  0,5887  +  0,33  (q  (0)  -  0,87} ,  (42) 

{0,87  < q  (0)  0,93  Bee.  flojnr], 

ca  (*„.)  =  0,6085  +  0,1758  (q  (0)  -  0,93}  -  4,29  (q  (0)  -0,93}*  (42a) 
[0,93 ^  q (0)< 0,99  eec.  flo/m], 

L  =  2,0+21,7  (q(0)  — 0,87}  +  3,88-10*  (q  (0)— 0,87}*,  (43) 

[0,87  ^q(O)^  0,93  Bee.  fiona]. 


L  =  4,7  +  24,17  {q  (0)  _  0,93}  +  2,04  •  10’  {q>)_0,9l}»,  (43) 

[0,93  ^  q  (0)  ^  0,99  boc.  ao:iii]. 

*==100,  [0,87  <^q  (0)^0, 93  nee.  aojhiJ.  (44  j 

*  -  ICO-  0,391  - 10* [q  (0)  -  0,93}  -  0,417 {q  (0>-0,93}*,  (44a) 
[0,93  ^q(0)<  0,99  BCC.  flO.-IItJ. 


i!,Note  to  formula  (39)*  ’lithin  the  brackets  are  indicated  the 
ranges  of  change  of  inputs  for  which  these  dependences  hold. 


The^  above  formulas  for  all  outputs  as  functions  of  the 
inputs  c, (0)  and  1(0)  show  that  the  plasma-chemical  reactor 
is  a  non-linear  transformer  for  these  inputs,  and  the  curves 
we  have  constructed,  together  with  the  corresponding  approxi¬ 
mate  analytical  expressions,  are  its  non-linear  characteristics „ 
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2.  The  Influence  of  fluctuating  initial  conditions  on 
the  solution  to  equation  set  Ho ) :  setting:  u'p 'the  problem. 


In  carrying  out  the  calculations  above,  the  initial  con¬ 
ditions  'Z.  (0)  and  T(0)  were  assumed  to  be  independent  of  time. 
In  connection  with  this  the  characteristics  obtained  for  the 
reactor,  when  the  reactor  is  regarded  as  a  non-linear  trans¬ 
former,  are  the  reactor's  static  characteristics.  Under  real 
conditions  the  initial  concentration  of  methane  and  the  ini¬ 
tial  temperature  of  the  gas  mixture  undergo  random  fluctuations 
over  time.  We  shall  describe  in  terms  of  stationary  random 
functions  (source  D0iD$Q)  these  mathematically  random  time 
fluctuations  of  the  initial  concentration  *0,  (0)  and  of  the  ini¬ 
tial  temperature  T(0),  assuming  (0)  anu  T(0)  to  be  absolutely 
continuous,  normal,  random  functions  of  . , *  Since  the  random 


*In  order  to  compare  the  results  obtained  for  such  an  approach, 
with  experimental  data,  it  is  evident  that  it  must  also  be 
assumed  that  these  random  functions  are  ergodic.  For  this  it 
is  sufficient  that  their  self-correlative  functions  do  not 
approach  zero  too  slowly  (source  Q3>. 


functiOTis  under  examination  vary  over  a  finite  interval 
Lx».  x*J  (where  x  is  either  of  the  two  functions  *5^(0)  and 
T(0;  it  is  necessary  to  use  truncated  normal  distributions 
(source  M>. 

The  one-dimensional  truncated,  normal  law  of  distribution 
has  the  following  form: 

4**  *!<*.. 

4  r  (45) 


where  x,  =  x(tj  )  and  6^  is  the  mean  square  deviation  of  t*,a 
random  rune t ion.  The  magnitude  £  is  chosen  such  that  f.  (x.) 
is  normal  to  unity  in. the  interval  of  definition  for  the  ran¬ 
dom  function. 

If  the  inputs  of  a  non-linear  transformer  are  stationary 
random  functions  of  time,  then  the  transformer's  outputs  Jj, 
'S'WzJ,  and  a  will  also  be  stationary  random  functions  of 
time  (as Burning  that  the  non-linear  transformer  has  stationary 
characteristics).  For  the  determination  of  the  statistical 
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characteristics  of  these  random  functions  we  use  formulas 
( 39 ) — (4^) .  Here  it  is  assumed  that  the  non-linear  transfor — 
mer  is  non- inertial  (sources  (321 »  £>a>;  that  is,  any 

output  y(t)  at  a  given  moment  or  time  t_  is  expressed  as  a  func¬ 
tion  of  any  input  x(t)  at  the  same  moment  of  time: 

»  W  “  f  ix  ($.  (46) 

3.  The  determination  of  the  mathematical  expectations 
and  of  the  correlation  functions  of  L7  w:  and  sr. 

The  basic  practical  consideration  in  the  study  and  con¬ 
struction  of  non-linear  systems  at  the  present  time  is  know¬ 
ledge  of  the  correlation  function  and  the  mathematical  expec¬ 
tation  of  the  random  output  function.  The  mathematical  expec¬ 
tation  <y>  of  a  stationary  random  process  y(t)  for  an  output 
of  a  non-linear  transformer  is  determined  in  the  following 
way  (source  M): 


<»>  =  <yi>  +  <yt>. 


(47) 


where 


<y»>  -  {2 o'  «*>  - x.)1-1  \Qt  {£,)-  <?!(£.)]>,} , 


Km  M 
*  ,  * 


(48) 


and 


k«0  Z«0 


)  *“»(— rf}* 


5  = 


*  -  <*> 


(49) 

'(50) 

(51) 


Here  A.  ,  are  the  constant  coefficients  in  formulas  (^9)- 
(44) ,  x£  is  the  value  of  £  dividing  the  region  of  variation 
,  x*]  into  two  equal  parts. 


The  computational 
of  the  output  n¥ ( e  ) 

Da ) ; 


i 


formula  for  the  correlation  function 
written  in  the  following  form  (source 
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(52) 


*v<e)-S  #«*>.<**)  4r-- 

uitat-e  "”l 

«.«*>.  3,)  -  B  S  S  ^1,  «i  «*>  -  *.)“[. «■„  (t.)  - 

+  S  ''*■  ilTTriiT  x 

>-e  i-t 

X  4«*>  -  x6)t  l  (A f,„  (VJ  -  Af,n  (go)], 

e  =  ti — ij,  (53) 
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where  rv  ( &  )  is  the  normal  self-correlation  function  of  input 

s 

^  dl,  (54) 

— • <0  * 

♦  (55) 


In  practice  in  computations  using  formuU  (52)  it  is 
necessary  to  be  limited  by  the  finite  number  of  terms  of 
order  m.  The  evaluation  of  accuracy  of  the  approximated  deter¬ 
mination  of  the  correlation  function  Xy(e)  is  given  by  the 
following  formula  (source  ) : 


\RV(Q)-Km<<z, 


1  an  i 

(56) 

i  r";e) , 

ft! 

(57) 

»*(<*>.  «,) 

n! 

(58) 

The  results  of  computing  the  dependences  of  the  mathe-  ' 
matical  expectations  ^  (z*^  ,  ('i^and  L  on  the  mean- 
square  deviation  of  initial  temperature  (for  a  given  mathe¬ 
matical  expectation  of  the  initial  temperature  \T{0)y  =  3000°K 
and  a  constant  value  of  the  initial  concentration  of  methane 
'SJ  (0)  =  0.95  by  weight)  are  shown  in  Figure  15a,  curves  1-3 
respectively.  Figure  15b,  curves  1-3  show  the  dependencies  of 
('Sj(Zrt1)/  ,  ,  and  {  L/ on  the  mean-square  deviation  of  the 

initial  concentration  of  methane  (where  the  mathematical  ex¬ 
pectation  of  tne  initial  concentration  of  methane^?!’,  (0)^  =  0.95 
by  weight  and  the  constant  of  the  initial  temperature 
T(o)  =  3500°;:  fsicn  ). 


<£>,  Cm  (Cj(7m)},  Dec  done  (J> 


Figures  15a,  15b. 
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For  Figures  15a,  15b:  The  dependence  of  the  mathematical 
expectation  of  (1)  the  maximum  concentration  of  acetylene 

,  (2)  the  degree  of  decomposition  of  methane (  , 

ancT  (3)  the  length  of  the  reaction  chamber  ^on  the  mean- 
square  deviation  of  different  inputs  to  a  plasma-chemical 
reactor  v(0)  =  3*10^  m/sec.  . 

Key:  2.  %  (zm) ,  by  weight 

1.  ^(z^),  by  weight  3.  %  by  weight. 

On  the  basis  of  these  curves  the  following  approximate 
analytical  functions  were  set  up: 

<e, (*J>  *=  0,475  —  3,166 •  10~*cjt  — 1,833.  iO"* a}, 

<L>  =  21  —  1 ,135  ■  10-*  aT  +  0,355  - lCr4  o’., 

<*>  =  83  —  0,217  •  10-*  oT  —  0,283 ■  i0~*  4 
0<or<300°K). 

<«*  (*.«)  >  =  0,6103  —  0,425  <jm  -|-  9,5  a],  , 

<£>  =  6,0+ 15<jm +5.10*3^, 

<»)  =  99,2  -40  a^  +  2  -10’  a£ 

0  ^  0,02  % 

A  calculation  of  formulas  (52),  (53) »  and  (5*0  gave  the 
following  formula  for  the  correlation  function  of  the  output 
&  (z«)  Ha  (0).  (In  the  calculation,  m  •  3,  T(0)  =  3000°K, 
c, (0)  »  0.95,  and  <£,.  =  100°K) : 

^(0)  =  1,122  lO-arT(0)  + 1,62- 10-*4(9> +6,61  10-»^  (0).  (61) 

The  relative  error  in  the  determination  of  the  correlation 
function  of  the  output  ^(z*,)  is  about  10$. 

4.  The  influence  of  the  values  of  the  constants  of 
reaction  rates  on  the  inputs  of  a  plasma-chemical  reactor. 
for  the  hinh-temperature  decomposition  of  methane  and  ethy¬ 
lene. 


(59) 


(60) 


It  is  viell  known  that  the  values  found  by  various  authors 
for  the  energy  of  activation  and  the  pre-exponential  coeffi¬ 
cients  of  the  constants  of  reaction  rates,  for  one  or  another 
chemical  reactions,  often  fluctuate  over  wide  ranges;  hence, 
selection  of  their  most  probable  values  is  a  matter  of  consi¬ 
derable  difficulty.  Such  a  situation  is  true  for  the  process 
of  conversion  of  methane  into  acetylene  (sources  I_44-46[3  ). 

Table  6  shows  the  reaction  process  and  values  for  the 
energy  of  activation  and  for  the  pre-exponential  coefficients 
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of  the  consta-  of  reaction  rates  which  have  been  obtained 
by  various  authors  for  this  process  in  different  experiments 
using  various  methods. 


PfJKUKH® 

**,  «i 

Ej  fC*T/I4/.H^>  | 

© 

JlHrrpiryp.i 

2CH.ic,H,  +  2H, 

4,5.10*» 

91 

|44) 

10“ 

103 

I45i 

Cifli*-*  CtHi  -f 

2,57-ln* 

40 

I«1 

it4 

2,2-10“ 

60 

146] 

CfHr-*  2G  -f  Hj 

1,7-10* 

30 

144) 

"table  6.  The  values  of  the  constants  of  reaction  rates 
of  the  process  of  the  pyrolysis  of  methane,  according  to 
various  sources. 

key: 

1.  reaction  3-  F  in  kilocal. /mole 

2.  k°  in  sec.-1  4.  "source. 

In  order  to  clarify  the  question  of  the  sensitivity  of 
the  mathematical  model  of  the  plasma-chemical  reactor  to  a 
change  in  the  constants  of  reaction  rates  k,  and  k„,  for  dif¬ 
ferent  initial  temperatures  T(0),  we  have  constructed  the 
static  characteristics  reflecting  the  output  $3  (zm)  as  a 
function  of  input  T(0)  for  various  selections  of.  the  indi¬ 
cated  constants  and  for  fixed  inputs  v(Q)  =  3*104  cm. /sec. 
and  (0)  =  0.9  by  weight.  For  this  purpose  an  integration 
by  computer  was  carried  out  on  equation  set  (36)  for  the 
following  combinations  of  constants  taken  from  Tabic  6; 

(kf,  Up,  Ckf,  kj),  (kf,  kp,  Ckf,  kp« 

and  for  the  following  three  inputs  of  T(0)  :  2.5  10-^; 

3.0*103  and  3.5-loJ,  for  fixed  inputs  v(0)  =  3*10^  cm. /sec. 
and  o(0)  =  0.9  by  weight. 

*Thn  indices  A,  2,  and  C  for  the  constants  of  react  ion"  "rates 
of  the  chemical  reactions  indicate  that  they  were  taken  from 
sources  D’-G,  C'-'53,  and  C  463  respectively.  The  constant  k« 
was  in  all  cases  the  same  (see  Table  C) . 


The  computations  were  carried  out  at  the  computer  center 
of  the  Academy  of  Sciences  of  the  I’SST  on  a  "Strela"  Arrow 
computer  using  the  Adams -S termer  method  vrith  automatic;, selection 
interval  of  integration  and  a  relative  accuracy  of  10”  ' .  The 
computation  time  for  one  variant  was  about  10  minuter,.  The 
statistical  characteristics  obtained  as  a  result  of  these  com- 
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putations  arc  shown  in  Figures  16a,  16b.  An  examination  of 
the  curves  shown  in  these  figures  permits  the  following  quali¬ 
tative  conclusions  to  be  made: 

1.  As  could  have  been  expected,  the  influence  of  a 
change  in  the  constants  of  reaction  rates  was  different  for 


Figures  16a,  16b.  The  dependence  of  various  outputs -of 
a  plasma-chemical  reactor  on  the  initial  temperature  T(0)  for 
different  selections  of  constants  of  reaction  rates  kp  and  k^,. 
16a:  the  maximum  concentration  of  acetylene,  %<*■>■  <  curves  I ) 
and  the  degrees  of  decomposition  of  methane,  §  (curves  II). 

1 6b :  the  length  of  the  reaction  chamber. 

Key : 

1.  shares  by  weight. 

2.  Replacing  the  constant  of  the  reaction  rate  of  the 
decomposition  of  methane  kf  with  kf  has  practically  no  effect 
on  the  values  and  intervals  of  the  curves,  but  replacing  the 
constants  of  reaction  rates  of  the  decomposition  of  ethylene 
(k*  for  k£  )  has  a  substantial  effect  on  the  curves.  This  is 
naturally  explained  by  the  fact  that  in  the  given  temperature 
range,  the  relative  change  k  is  small  compared  to  the  relative 
change  of  k^. 

3.  In  accordance  with  the  physical  sense  of  the  constant 
of  reaction  rate  of  a  chemical  reaction,  a  change  in  k^_  most 
of  all  affects  the  effective  length  of  the  reaction  chamber 
(that  is,  in  practical  terms,  it  affects  the  time  of  the  reac- 


tion  up  to  the  point  when  the  formation  of  reaches  a  maxi¬ 
mum).  It  affects  less  the  absolute  magnitude  of  this  maxi¬ 
mum.  Hence  it  follows  that  the  experimental  value  of  kji  must 
be  determined  with  far  greater  accuracy,  and  the  possible 
sources  of  errors  and  discrepancies  in  the  values  of  ob¬ 
tained  in  sources  CMO  and  {f  463  should  be  analyzed. 

In  order  to  clarify  the  effect  on  the  numerical  solution 
of  equal  relative  changes  of  k.  and  k»,  we  examined  the  inter¬ 
mediate  value  of  the  latter  constant  k-  --  5-8*10°  exp(-5l/RT), 


for  which 


*? 
*?  ’ 


that  is,  its  relative  change  is  equal  to  the  relative  change 
of  k,  . 

Figures  l6a  and  16b  show  the  results  of  the  computation, 
given  equal  relative  changes  in  k(  and  k^.  As  could  have  been 
expected,  the  effeeb  of  replacing  k?  with  k*  is  somewhat  larger 
than  that  of  replacing  kf  with  kf  (the  relative  change  of  the 
constants  in  both  cases  being  equal).  The  influence  on  the 
time  of  the  reaction  is  not  large  and  decreases  with  an  increase 
in  temperature.  The  influence  on^(zm)  is  somewhat  greater 
for  higher  temperatures  than  for  lower  ones,  and  greater  than 
the  influence  of  an  equal  relative  change  of  kf  for  kf  . 

It  is  possible  that  the  results  found  above  have  a  some¬ 
what  more  general  significance  for  the  high-temperature  processes 
that  are  a  successive  link  in  the  chain  of  conversions.  In  each 
concrete  case  the  decisive  role  will  of  course  be  played  by  the 
ratios  of  the  magnitudes  of  the  constants  and  the  "remoteness" 
of  the  reaction  to  which  they  relate  from  the  end  product 
whose  output  is  being  monitored.  At  the  present  time  there 
is  no  criterion  to  indicate  to  what  degree  this  or  that  constant 
of  reaction  rate  affects  the  output  of  the  plasma-chemical 
reactor.  The  results  obtained  point  to  a  possibility  of  using 
the  following  sensitivity  criterion  of  the  outputs  of  the  plasma- 
chemical  reactor  to  a  change  in  a  constant  of  reaction  rate: 


l£*  (*Tn)IkX  [£»  (zm)])cii 


M£m)Jn 


[‘-Srf. 


where  (z^)  J  kx ,  £  c4  (z«n)  k-  are  the  values  of  output 
cs  ( Zff, )  corresponding  to  two  different  values  of  the  same  con¬ 
stant  of  reaction  rate,  kx  and  k^  respectively. 

It  is  possible  to  construct  analogous  criteria  of  sensi¬ 
tivity  for  the  outputs  L  and  £, 


CO 
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Tn  conclusion  we  note  that  the  question  touched  on  above 
is  a  part  of  a  ..ore  general  problem  concerning  the  sensitivity 
of  solutions  of  an  equation  set  of  chemical  kinetics  to  the 
change  of  the  constants  of  reaction  rates  that  are  embodied  in 
the  equation  set  (see  Section  V,  beginning) . 

The  investigation  that  has  been  carried  out  allowed  us, 
to  construct  the  static  characteristics  of  the  process  and  to 
show  that  the  reactor  may  be  regarded  as  a  non-linear  trans¬ 
former.  In  the  latter  sense  approximate  analytical  dependences 
may  be  established  between  the  inputs  (for  instance,  the  ini¬ 
tial  concentration  of  methane  (0),  the  initial  temperature 
?(0)  of  the  gas  mixture)  and  the  outputs,  that  is  to  say,  the 
magnitudes  characterizing  the  process  of  conversion  of  methane 
to  acetylene  (for  instance,  the  maximum  concentration  of  the 
end  product  ^(z^),  the  effective  length  of  the  reaction  cham¬ 
ber  L,  and  the  degree  of  decomposition  of  methane,  _§),  These 
dependences  may  be  established  for  several  selections  of  con¬ 
stants. 

The  relationships  found  for  reactor  outputs  as  functions 
of  reactor  inputs  enable  the  clarification  of  the  degree  of 
sensitivity  of  various  outputs  to  one  particular  input,  and 
also  the  effect  of  one  particular  input  on  various  outputs, 
both  for  the  case  of  fixed  inputs  as  well  as  the  case  of  con¬ 
tinuous  random  fluctuations  of  the  inputs.  The  results  ob¬ 
tained  for  this  enable  a  comparison  of  the  effect  on  the  out¬ 
puts  of  fluctuating  inputs  with  the  effect  on  outputs  of  inac¬ 
curacy  in  the  given  constants  of  reaction  rates  used  in  the 
equations.  In  particular,  for  the  process  of  the  pyrolysis  of 
methane  examined  above,  it  can  be  seen  from  a  comparison  -of 
curves  i  in  Figures  14a  and  14b  that  in  the  given  region  of 
change  for  the  magnitudes  of  inputs  T(0)  and  (0),  the  maximum 
of  the  concentration  of  acetylene  (Zfn)  is  more  sensitive  to 
a  change  in  input  T(0);  it  is  evident  that  this  must  be  taken 
into  account  in  the  regulation  of  output  'hi  (z*),  In  the  case 
where  the  inputs  f  in  our  case,  T(Q)  and  '3|  (o)j  continuously 
and  staticnarily  fluctuate,  they  may  be  mathematically  modeled 
with  the  aid  of  stationary  random  functions  of  time.  The  ap¬ 
proximate  static  characteristics  of  the  process  enable  us  to 
find  the  mathematical  expectations  of  the  outputs  (in  this 
example,  ^(z^),  L  and  s,}  as  functions  of  the  mean-square  devia¬ 
tions  o'’  the  inputs — in  this  example,  T(0)  and  o'.  (0) — and  also 
compute  the  self-correlation  functions  of  the  outputs — for 
example,  ^(z^).  In  particular,  the  graphs  of  mathematical 
expectaf  on  <  ^  (zm)^  in  Figures  lya  and  15b  show  that  for  ran¬ 
dom  oscillations  of  inputs  (0)  a;F  T(0),  the  average  value 
of  the  output  of  the  end  product,  acetylene,  “03(2^)  may  be  as 
much  as  about  5-5u  lower  than  the  value  of  this  magnitude  for 
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constant  initial  conditions. 

The  static  characteristics  of  the  plasma-chemical  reactor 
that  have  been  obtained  in  this  article  may  be  used  for  the 
solution  of  a  number  of  problems: 

a)  for  investigating  the  question  of  controlling  the 
end-condition  of  the  processCe.g.  achieving  a  maximum  of  the 
concentration  of  acetylene  by  means  of  changes  in  the  inputs 
*&l(  0)  and  T(o)l; 

b)  for  solving  the  more  general  problem  of  choosing  the 
optimal  (for  a  given  criterion  of  optimality)  value  of  the 
inputs — e.g. ,T(o)  and  0^(0)  and  others. 

The  correlation  functions  of  outputs  and  inputs  may  be 
used  for  solving  problems  of  predicting  the  behavior  of  out¬ 
puts  of  the  reactor  given  fluctuating  inputs  {by  using  the 
theory  of  prediction  of  stationary  random  functions). 


It  can  be  seen  in  Figure  16a  that  the-  scatter  of  values 
for  the  output  £i(z  )  caused  by  inaccuracy  in  the  determi¬ 
nation  of  constantsmof  reaction  rates  k-,  and  kp,  for  different 
values  of  input  T(0),  lies  in  the  2-10$  range.  ■  Since  the 
decrease  in  output  'cd(z  )  due  to  random  fluctuations  of  in¬ 
puts  is  of  about  theJsaffle  order  of  magnitude,  these  fluctu¬ 
ations  must  be  taken  into  account  in  any  calculation  for  the 
plasma-chemical  reactor. 


The  approach  developed  in  this  article  may  be  applied 
without  difficulty  to  problems  with  a  great  number  of  in¬ 
puts,  outputs,  and  constants  of  reaction  rates'. 

It  should  be  noted  that  we  have  examined,  strictly  speak¬ 
ing,  one  section  of  the  characteristic  surface  in  that  multi¬ 
dimensional  space  called  "input-output;"  to  construct  approx¬ 
imately  that  surface  in  its  entirety  could  be  done  only  with 
the  expenditure  of  a  large  amount  of  computer  time. 

In  conclusion  we  note  that  certain  of  the  problems 
examined  above  are  problems  of  information  theory  as  well, 
and  warrant  the  corresponding  written  exposition  in  that 
field  as  well  . 


V.  The  influence  of  the  values  of  the  constants  of 
reaction  rates  of  chemical  reactions  on  the  solution  of  a 
set  of  equations  of  chemical  kinetics  (for  example,  the  oxi¬ 
dation  of  methane).  "  ' 

In  looking  at  one  or  another  variant  of  a  reaction 
mechanism,  chemists  evaluate  the  possible  influence  of  a 
change  in  value  of  certain  constants  of  reaction  rates  on 
the  course  and  results  of  the  reaction,  .if  the  reactions  are 
relatively  non-complex.  It  is  evident  that  such  evaluations 
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cf  the  sensitivity  of  a  reaction  mechanism  to  a  change  in 
values  of  constants  of  reaction  rates,  if  such  evaluations  are 
even  possible,  become  more  difficult  in  the  case  where  quite 
complex  reactions  are  under  study.  It  is  in  this  regard  that 
the  problem  arises  of  working  out  objective,  quantitative  cri¬ 
teria  of  "sensitivity"  of  mechanisms  of  complex  chemical  reac¬ 
tions  and  standard  methods  for  computing  these  criteria. 

As  will  be  shown,  the  method  developed  earlier  (see  Sec¬ 
tion  I)  allows  us  to  obtain  by  computer  useful  and  objective 
information  about  the  sensitivity  of  mechanisms  under  study 
to  changes  in  different  constants  of  reaction  rates  in  various 
stages  of  the  reaction,  for  various  conditions  under  which  the 
reaction  is  proceeding. 

This  section  is  devoted  to  the  working  out  of  a  mathema¬ 
tical  method  which  will  permit: 

1)  giving  a  quantitative  expression  of  sensitivity  of 
this  or  that  reaction  mechanism  (for  a  given  aggregate  of 
constants)  to  a  change  in  values  of  the  constants  of  reaction 
rates  of  individual  stages,  using  the  quantitative  criteria 
of  sensitivity  to  be  introduced  below: 

2)  constructing  an  algorithm  of  computation  of  the  cri¬ 
teria  we  have  introduced; 

3)  finding  the  dependence  on  time  of  these  criteria  for 
various  conditions:  under  which  the  reaction  is  proceeding  and 
for  various  selections  of  constants. 

This  method  will  aid  in  increasing  the  information  derived 
from  experimental  kinetics  data,  and  will  help  give  the  chemist 
objective  criteria  for  evaluating  and  comparing  the  influence 
of  various  constants  in  various  stages  of  this  or  that  reac¬ 
tion. 

The  methodology  proposed  is  illustrated  below  for  the 
example  of  the  oxidation  of  methane. 

1.  Formulation  of  the  nnnhl em. 

As  is  well  known  the  problems  of  chemical  kinetics,  from 
a  mathematical  point  of  view,  come  down  to  the  problem  of  sol¬ 
ving  and  analyzing  sets  cf  non-linear,  ordinary  differential 
equations  of  the  first  order  with  coefficients  representing  • 
constants  o.f  reaction  rates  k~(q  =  l,2,..,,m).  For  the  sake 
of  simplicity  the  exposition  That  follows  examines  the  case 
of  «  co thermal  reactions  (for  which  all  ka  =  const.)  under 
"to tie  conditions;  expanding  the  results  "obtained  below  to 
non- isothermal  reactions  in  a  current  stream  does  not  offer 
-1 '  ff! unities  in  principle. 

The  cf  equations 


of  chemical  kinetics  may  be  written 


in  the  following  form: 

'HL  =  /i(c-  *)(/  =  *,  2 . n)  (62) 

for  initial  conditions  0^(0)  »  c9,  where  the  function 
fj(o,k)  is  the  sum  of  terms  of  the  form 

■‘qfrlPM’. 

where  c.=  arc  the  present  concentrations;  n  is  the  number 
of  components  taking  part  in  the  reaction;  jo,s  is  the  order 
of  the  reactions. 

The  solutions  to  equation  set  (62)  are  functions  of  time 
and  depend  on  k.  as  well  as  on  the  parameters,  namely- 
c .  =  c.(t,k).  From  the  linearity  of  the  right-hand  sides  of 
equation  set  (62)  relative  to  kq,  it  follows  (source  l47J) 
that  the  derivatives  3  c.(t,k)/a  kQ  exist  and  are  continuous. 
Hence  the  mixed  second-order  derivative  of  the  functions 
c,(t,k)  with  respect  to  t  and  k  exist  and  are  continuous. 
Differentiating  the  equation.  se$  (62)  with  respect  to  kq 
(for  a  given  aggregate  of  kq)  and  designating  H 

uh(^k)  —  • 


we  obtain  the  following  set  of  equations  for  the  functions 


gif  = 

whe**^ 

**(«(<),  «(/).*!  = 

dki 

k)  -  JlteML  , 

ocl 


(63) 


n 


+  2  anlc(lh  kj  U[„  [eh),  A], 


In  view  of  the  fact  that  the  functions  F;q  depend  not 
only  on  Ujq  but  also  directly  on  c.,  the  equation  sets  (62) 
and  (63)  should  be  solved  concurrently  for  the  following  ini¬ 
tial  conditions: 

* 

c/  (0)  =  c®  ;  (0)  =  0.  (64) 

The  second  of  the  two  equalities  in  conditions  (64)  reflects 
the  fact  that  c®  does  not  depend  on  the  constants  of  reaction 
rates  It  .  J 

T. 

Solving  concurrently  the  equation  sets  (62)  and  (63), 
we  obtain  all  the  concentrations  c.  and  also  functions  u. 
for  all  moments  of  time.  At  present,  it  is  not  possibleJclto 
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solve  such  equation  sets  (62)  and  (63),  even  for  not  very  com¬ 
plex  chemical  reactions,  except  on  electronic  computers.  The 
method  proposed  in  section  I  of  this  article  enables  us  to 
obtain  a  numerical  solution  to  the  equation  sets  on  a  compu¬ 
ter,  to  the  required  degree  of  accuracy. 

The  basic  difference  between  equation  set  (63)  and  equation 
set  (62)  is  that  the  functions  u^_  may  take  both  positive  and 
negative  values  (in  distinction  from  c.),  and  also  that  there 
are  no  limits  whatsoever  on  the  magnitudes  of  Uja.  These 
peculiarities  of  equation  set  (63)  are  not  crucial  if  the 
numerical  integration  is  carried  out  with  the  use  of  the  algo¬ 
rithm  worked  out  earlier  in  this  article. 


However,  it  should  be  taken  into  account  that  for  the 

ujq  as  well  as  for  the  concentrations  e.,  the  equa¬ 
tion  of  balance  must  be  satisfied.  Actually  at  any  moment 
of  time  t  Huiq  (t)  =  0,  and  inasmuch  as  HE u*a(0)  »-0  from 
conditions  (64;,  u*a(t;  =  0.  Thus  solving  one  set  of  equa¬ 
tions  (62)  and  (63),Jin  addition  to  the  information  derived 
from  the  solution  of  the  equation,  set  (62)  of  .  chemical  ki¬ 
netics,  also  yields  quantitative  evaluations  of  the  sensiti¬ 
vity  of  the  solutions  of  equation  set  (62)  to  a  change  in 
the  constants  of  reaction  rates  embodied  in  the  equation  set. 


2.  Criteria  of  the  sensitivity  of  equations  of  chemical 
kinetics  of  isothermal -reactions  to  a  change  in  the  constants 
of  reaction  rates  embodied  in  the  equation  sek. 

Using  the  numerical  solution  of  equation  set  (62)  and 
(63)  it  is  possible  to  determine  the  following  criteria  of 
sensitivity  of  the  set  of  equations  of  chemical  kinetics  (62) 
to  a  change  in  the  constants  of  reaction  rates  embodied  in 
it. 


Most  evident,  probably,  but  not  the  most  suitable  cri¬ 
terion  is  the  very  function  itself: 

(65) 

The  difficulty  in  using  this  criterion  is  connected  with  the 
known  fact  that  kq  may  take  on  various  dimensionalities 
because  of  the  various  orders  of  the  reactions  that  are  a  part 
of  the  mechanism  under  study.  Also  the  values  of  k  may  dif¬ 
fer  from  each  other  by  many  orders  of  magnitude.  Moreover, 
the  various  concentrations  c^  may  differ  greatly  from  each 
other,  as  vrell  as  change  considerably  during  the  course  of 
the  reaction.  Because  of  all  of  the  foregoing  it  is  expe¬ 
dient  to  introduce  the  following  criteria  of  sensitivity:* 


(66) 


W '■*>  _  ku  a  ,, 
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(67) 

(68) 
(69) 


■"'ve  note  that  criteria  (6?)  and  (68)  may  be  used  only” when 
values  of  ci  4  0. 

tj 


where  cVlc)  is  some  characteristic  (depending  on  the  con¬ 
stants  of  reaction  rates  k  as  well  as  on  the  parameters) 
value  of  the  concent  rat  ion^-of  the  i-th  component.  The  maxi¬ 
mum  value  of  c*  may  be  taken  as  #■?!<). 

J  J 

The  criteria  (65) -(69)  give  answers  to  several  different 
questions.  For  example,  criterion  (65)  shows  (for  all  values 
of  t,)  to  what  extent  the  concentrations  c .  change  as  a  result 
of  a  change  of  A  k  in  a  constant  of  reaction  rate  k  ,  where 
6k  — »  0.  At  the  Same  time  criterion  (68)  determines  for 
each  moment  of  Lime  the  relative  changed c-/c.  in  the  con¬ 
centration  Ci  caused  by  the  corresponding  change  &  kg/ky—XD 
of  the  constant  k  ,  etc.  The  constants  of  sensitivity  C65)- 
(69)  rnay  be  particularly  useful  in  evaluating  the  influence 
of  a  change  in  each  of  the  constants  k  on  any  of  the  con¬ 
centrations  c.  for  any  moment  of  time  c,  and  also  for  evalu¬ 
ating  the  infiuence  of  any  of  the  constants  on  all  or  -some 
of  the  concentrations  for  any  moment  of  time.  It  is  evident 
that  in  the  course  of  this  procedure  it  becomes  clear  which 
of  the  constants,  and  for  what  stages  of  the  reactions,  (and 
at  what  moment  in  time)  are  crucial  or  conversely,  have  lit¬ 
tle  influence  on  the  reaction  (for  a  given  aggregate  of  k  ). 
Such  research  may  be  carried  out  for  any  initial  ratios  or 
the  concentrations  and  any  initial  thermodynamic  conditions. 

17e  shall. make  one  more. remark.  All  of  the  results  which 
have  been  described  concerning  the  criteria  of  sensitivity 
are  good  for  infinitesimally  small  changes  in  the  constants 
kg  (for  a  given  aggregate  of  k  ) .  In  practice,  as  it  is  usu¬ 
ally  done,  the  final  increments  of  kq  are  examined.  Inas¬ 
much  as  the  solutions  of  equation  set  (62)  and  (63),  yiel¬ 
ding  Cj  =  Cj(t,k)  and  =  Uj„(t,k)  are  for  most  moments 
of  time  smooth  functions^of  the  parameters  k  ,  such  an  assump¬ 
tion  does  not  lead  to  significant  errors  if^the  necessary 
v/atchfulness  is  maintained.  It  is  evident  that  a  sufficient 
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mathematical  condition  for  this  is  the  smallness  of  higher- 
than-first-order  terms  in  the  expansion  of  the  functions 
c<(t,k)  by  degrees  6ka  in  size.  Since  equation  set  (63; 
is  linear  with  respect  to  u<q,  the  sensitivity  of  the  con¬ 
centrations  to  a  simultaneous  change  of  several  constants  kq 
may  be  shovm  as  the  sum  of  the  respective  sensitivities  of 
the  concentrations  c<  to  a  change  of  the  individual  constants 
£taken  separately^. 

3.  The  reaction  of  the  oxidation  of  methane^. 

As  an  example  of  the  use  of  the  method  outlined  above  we 
shall  examine  the  isothermal  reaction  of  the  gaseous  oxidation 
of  methane.  The  reaction  process  of  this  reaction,  not  inclu¬ 
ding  the  reactions  of  chain  breakage,  is  described  in  sources 
OSjf  and  £493  and  is  shown  in  Table  7.  In  Table  7  thevaiues 
of  the  constants  ka  for  T  =  103  °K  are  converted  to  the  dimen¬ 
sion  sec.-*  using  the  following  relationships: 

kt  =s  ktc0\  c.  =  — =  0,122 -10'*  m  oL t/c'm1 

^  /  0»  d 

T0  =  273°  K .  V0  =  2,24-10*  cafymde. 


{ 

s. 

<S> 

I  : 

E«. 

3na<UBnn  7/ 

(T-1000  *K) 

1 

CH.  +  O^CHj  +  HO, 

Iff* 

55 

0,H'610-« 

CH,-f-Oj-.  CH,00 

low 

0 

0, 122.10’ 

CH»00—  CH,0  4-  OH 

10”** 

20 

0,425-10* 

CH.  +  OH-  CH,+  H,0 

10M 

8.5 

0,168.10* 

CH,0  +  OH-  H.O  4-  HCO 

10” 

6,3 

0.512.1CX 

CH,0  +  Or-  HCO  4-  HO, 

10” 

32 

0,123.10* 

HCO  4-  0^  CO  4-  HO, 

10" 

0 

0,122-10’ 

CH,4-HOr-  H,0,  +  CH, 

10” 

19,75 

0,588-10* 

CH,04-H0r*  H,0,  4-  HCO 

10" 

8.5 

0,169-10* 

CO  +  OH-  CO,  4-  H 

10"  • 

7 

0,359- 10“ 

CH,  4-  H  -  CH*  4-  H, 

1.5-101* 

It, 2 

0,643-10" 

CH.0  4-H—  HC0  4-H, 

4-10"  Y~ 

2 

0,563-10* 

i  •  *,  -  otn  «*p  (-  e/rt).  ••  [»«i  - 


'T’ablc  7.  The  mechanism  of  the  oxidation  of  methane. 
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The  equation  set  (62)  for  this  reaction  and  the  part  of  equa¬ 
tion  set  (63)  corresponding  to  a  change  in  only  or*constant 
of  reaction  rate,  klf  has  the  following  form: 


~~  —  — ex  (Vj  4  Vt  4  Vi  4  Vci») 

—  —  Cj  (Vi  4  Vs  +  V«  +  V») 

**jp  —  —  Vi^s  +  cx  (Vs  4  Vs  4  Vc»  4  Vi) 
~~  =  —  (Vi  4  Vi)  4  C|  (Vi  4  Vi  4  Vi) 
•jp  =  —  Vi  4  V*£» 

-jp  =  —  ci  (V?  4  V*  4  Vi  4  ^ij^h)  4  Vi 
==  — Cj(Vi  4  V»  4  Vv)  4  Vi 
—  Ci  (Vi  4  Vi) 

pgr-  =  —  Vi c»  4  c*  (  Vj  4  V*  4  V«  4  kitCii) 

*jp  —  — kufiiclQ  4  V*c» 

^jp  =  c«  (Vi  4  Vi) 

~jp  =  C||  (VC1  4  Vc«) 

“jp  —  ^lO^lO 

-jp  ~  —  *uc,Cj|  —  kltCfCu  4  ^io^io 


(70) 


fir  (fir)  ~  ~  fir  ^‘C|  +  k*Ci  +  fc,c*  +  *»*»>  ~ 

~  c»  (fci  Hr + k*  fir  +  *•  fir  +  *“  Hr)  ~ ClC* 
fir  (fir) =  —  Hr(**c»  +  Av*  +  **Ci+*’c,)—  1  <71> 

~ c*  (*»  fir + *»  Hr + *•  it + fe?  Hr)  ~  c‘c* 
fi-  (Hr)  =  H7  +  ^*C<  +  *ucu  +  ktct)  + 


+  c»  VA‘»  5*r  d*7  ^  *“  SET 


4  -  (ft)  ■ - - fir  <*•“ + **•>  “  “  (A*  Hr + *•  5-) + 

.+  •gj7  (^ici  +  ^«e»  +  k-fit)  + 

+ c*  (*»  Hr + *■  "Sr +  *»  Hr) + c,Cj 

rf  /  dfj  \  ..  &•»  ,  L  -  •*»  ,  u.  dci 

nr\-m )  =  — A»  **r  + AiC*  w  +  A,Ca  **r 

»■  tk) -  ~  k  <** + **• + v< + - 
-«•(»•  fe+‘-fe.+*-t+*"Sr)+,:-w 

■3T  ( w)  ”  -  If  <**  + 1,',+  _ 

M&)  =  'fe(Vl  +  *,C,)  +  '’(*‘^  +l“&) 

f  (-fe)— **TK-**-fe-+  ' 


*v»ry  »i  •*»»  ■ 

-j-  -^-  (AjC7  -f  ^'«ci  +  + 

+ c»(*»  Hr  + A*  ■&  +  A*  "Sr +  *u  ft) 

fi(Hf)  — *«*  £ +**  4-  +**£ 
4r  (Hr) =  fir  +  Av,) + c<  (A*  fir + kt  fir) 

fi-(Hf) = Hr (A,iCi + ^ + Cm  (Aa  Hr + *1S  fir) 
fiKHf)^^’  Hf  +  A’,oC,#  fir 
fir  (Hr) =  ”  Hr (A,iC* + A,sCi)  ~ 

— «i4  (*u  +  *n  fir)  +  A‘»'<C;  ■sir  +  Ai#Cw  5fcT  j 


68  - 


Used  here  for  the  concentrations  (In  volumetric  shares)  of  the 
fourteen  components  taking  part  in  the  reaction  are  the  desig¬ 
nations  introduced  in  Section  III-l.  In  order  to  determine 
the  influence  of  the  magnitude  of  constant  on  the  solution 
of  equation  set  (70),  the  equation  set  (?Q)  and  (71)  must  be 
solved  concurrently.  In  the  case  where  it  is  necessary  to 
find  the  influence  of  a  change  in  other  constants,  equations 
analogous  to  those  of  equation  set  (71)  must  be  written  and 
then  solved  concurrently  with  equation  set  (70).  Thus  it  is 
possible  to  find  the  sensitivity  of  all  concentrations  with 
respect  to  all  constants  kq.  In  order  to  obtain  data  about 
the  sensitivity  of  all  concentrations  c .  to  a  change  in  one 
constant,  it  is  sufficient  to  write  an  Equation  set  of  the 
type  (71)  and  then  integrate  it  concurrently  with  equation 
set  (70).  However,  in  order  to  determine  the  sensitivity  of 
one  concentration  to  a  change  in  all  the  constants  k„  it  is 
necessary  to  obtain  a  solution  to  all  n(m  +  1)  equations  for 
functions  c *  and  u^q.  In  order  to  obtain  such  a  solution  it 
is  possibleJ(if  the  capacity  of  the  computer  permits  it)  to 
construct  and  integrate  concurrently  all  n(m  +  l)  equations. 

If  it  impossible  to  do  it  all  at  once,  however,  then  the  func¬ 
tions  u^q  =  b  Cj/?>  kq  should  be  integrated  first  for  one 
group  of  constants,  then  for  another,  and  so  on  until  all  the 
constants  kq  have  been  included. 

The  equation  set  (70)  and  (71)  for  various  k  was  inte¬ 
grated  on  the  computer  ‘BESM* at. the  computing  center  of  the 
Academy  of  Sciences  of  the  USSR  by  the  method  proposed  in  Sec¬ 
tion  1-3  of  this  article  with  a  relative  accuracy  of  from  10-^ 
to  10-5,  using  automatic  selection  of  the  interval  of  integra¬ 
tion.  The  computation  time  for  one  variant,  corresponding  to 
changing  one  of  the  constants  of  reaction  rates,  is  about  30 
minutes. 

In  the  integration  process  on  the  computer,  significant 
time  is  wasted  on  the  computation  of  those  concentrations 
L having  ci(0)  =  QJ  whose  values  become  greater  than  zero  very 
slowly.  ^Far  example,  for  an  integration  up  to  t  =  10~-  sec, 
about  50.5  of  the  machine  time  is  spent  on  the  integration  of 
equation  set  (70)  and  (71)  during  its  induction  period  (up  to 
t  a  10”'’'  sec.).  During  this  period  the  state  of  the  system 
has  practically  not  changed  at  all  from  that  of  the  initial 
conditions.  The  computation,  time  may  be  considerably  decreased 
by  means  of  removing  the  accuracy  check  with  respect  to  these 
concentrations  while  they  are  still  small  within  the  limits  cf 
accuracy.  In  the  computation  process,  balance  with  respect 
to  e<  as  well  as  tc  ujq  is  satisfied  to  an  accuracy  of  five 
or  six  significant  figures. 

A  discussion  of  the  results  of  the  integration  of 
equal ion  set  (70).  (71)  for  the  react  ion  of  the  oxidation  of 
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The  curves  shown  in  Figures  17-22  illustrate  the  compu¬ 
tations  for  criteria  of  sensitivity  ( 65)  —  ( 69 >  for  the  oxida¬ 
tion  of  methane.  These  curves  were  obtained  as  a  result  of 
the  integration  of  the  set  of  differential  equations  (?0), 
(71)  for  the  model  isothermal  reaction  of  the  gaseous  oxida¬ 
tion  of  methane  under  the  following  conditions: 

* 

T;  =  103  °K;  c,  (0)  ~  0,29^,  (0)  =  0,71;  Cp  (0)  =  0 

<p“3’4 . 14?:  ($■)<«,  5=0  .  14:  ?=*  1,2,...,  12),  (72) 

and  also  under  the  condition  of  a  fixed  aggregate  of  values 
of  constants  of  reaction  rates  for  the  reaction  (see  Table  7) 


Figure  17.  The  dependence  on  time  of  the  sensitivity 
of  the  concentration  of  formaldehyde  (05)  to  a  change  in 
constants  kq  (q  =  1,2,... ,12). 

Solid  lines:  dc^/  "bln  kq;  Dotted  lines:  "Vln  04/  *bln  kq. 

The  numbers  on  the  curves  correspond  to  the  numbers  of  the 
constants  kq  that  are  being  allowed  to  vary. 

The  line  of  mixed  dots  and  dashes:  the  concentration  of  for¬ 
maldehyde. 


Figure  18.  The  dependence  on  time  of  the  sensitivity 
of  the  concentration  of  water  (eg)  to  changes  in  the  con¬ 
stants  k4(  k6,  kg,  1<10. 

Solid  lines:  “&Cg/  "bln  kq;  Dotted  lines:  "bln  eg/  bln  kq. 

The  numbers  on  the  curves  correspond  to  the  numbers  of  the 
constants  k  that  are  being  allowed  to  vary. 

The  line  of  mixed  dots  and  dashes:  the  concentration  of  water. 
Key:  t,  in  sec. 
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Figure  19.  The  dependence  on  time  of  concentrations  of 
methane  (c^)  and  oxygen  ( C2 )  and  their  sensitivity  to  a  change 
in  the  constants  kg,  k^,  kg. 

Solid  lines:  ©c^/  ©In  k^;  Dotted  lines:  ©eg/  ©In  kq. 

The  numbers  on  the  curves  correspond  to  the  numbers  of  the 
constants  that  are  being  allowed  to  vary. 

The  line  of  mixed  dots  and  dashes:  ©In  c j/  ©In  kq. 

Key:  1.  t.  in  sec. 


Figure  20,  The  dependence  on  time  of  the  concentration 
S'-t,  (c-,1  and  its  sensitivity  to  a  change  in  constants  kzj,,  ks. 

-\q« 

9  olid  lines:  ©In  k^;  Dotted  lines:  ©In  c^/  ©In  kg. 

The  numbers  on  the  curves  correspond  to  the  numbers  of  the 


- 


constants  k^  that  are  being  allowed  to  vary. 

The  line  of  mixed  dots  and  dashes:  the  concentration  of  CH-j  (c^). 
Key:  3.  t  in  sec. 


Figure  21.  The  dependence  on  time  of  the  concentrations 
CH-,  (c3T,  CH-00  (ct),  and  HCO  (cq)  and  their  sensitivity  to  a 
change  in  the  constant  k^.  7 


Solid  lines:  ^c^/c^  l»ln  kj.;  Dotted  lines:  "ftc^/c^* In  k^. 

The  lines  of  mixed  dots  and  dashes:  the  concentrations  of 
CH3  °3  *  CH300  (c5)*  and  HC0  (09) • 

Key:  1.  t  : 


Figure  22.  The  dependence  on  time  of  the  sensitivity  of 
the  concentrations  Cl'  (c-J,  CH,00  (cj,  OH  (c 7),  HCC  (cQ) 
to  a  change  in  constant  kg.  J  y  ‘  y 
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NOT  REPRODUCIBLE 


The  numbers  on  the  curves  correspond  to  the  numbers  of  the 
concentrations  c^  under  study. 


Key:  1.  t  in  sec. 

In  Figures  17-22  are  shown  the  dependences  on  time  of  the 
criteria  of  sensitivity  of  five  concentrations  »  l°2l. 

(ClUl  ,  ICII20]  and  g20j  to  a  change  in  various  constants  or 
reaction  rates  (see  Figures  17-20);  as  well  as  the  criteria 
of  sensitivity  of  various  concentrations  to  a  change  in  two 
constants  of  reaction  rates  (k«  and  kg)  (see  Figures  21,22). 
For  several  of  these  figures,  ror  purposes  of  comparison,  are 
shown  curves  of  the  dependence  on  time  of  various  criteria 
of  sensitivity  for  the  same  concentrations  (cj  or  constants 
(k  ),  and  also — for  convenience— the  curves  of  the  dependence 
on-time  of  the  concentrations  themselves. 

As  can  be  seen  from  the  figures ,  the  values  of  the  cri¬ 
teria  of  sensitivity  may  be  positive,  negative,  or  approach 
zero.  The  positive  (negative)  value  of  a  criterion  of  sensi¬ 
tivity  at  the  moment  t  means  that  for  the  given  values  of  the 
constants,  an  increase  in  the  given  constant  of  reaction  rate 
leads  to  an  increase  (decrease)  of  the  value  of  concentration 
Ci  for  that  given  moment  of  time.  It  is  evident  that  when  the 
value  of  the  criterion  is  zero  at  time  Jt ,  this  means  that 
at  that  point  in  time  a  small  change  in  k  does  not  lead  to 
any  change  in  cj.  3 

An  examination  of  Figures  17-22  shows  also  the  following 
details. 


'--O  - 


rosj- 
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or  concentrations  having 
igur  :  ':'  17a,  l?b) ,  tas 
:c c t  to  a  change  in  various 


a  maximum 
curves 


for 


constant: 


’’l  into  two  categories:  a)  curves  that 


(for  example,  for 
the  sensitivity  with 
of  reaction  rates 

w  _  . .  _  _  remain  in  the  posi- 

.ivs  region  (or  the  negative  region)  of  values  of  the  criterion 
of  sensitivity  for  all  values  of  t  (that  is,  for  the  whole 
duration  of  the  reaction);  b)  curves  that  may  go  from  a  region 
of  one  sign  for  the  values  cf  tho  criterion  to  a  region  of  the 
ether  sign. 


^ov  concentrations  whose  values  simply  rise  continuously 
or>  fall  continuously  over  time  (for  example,  H«0  in  Figure  18), 
only  curves  of -the  first  category  occur.  * 


Tims,  the'  magnitude  and  s i«;n  of  the  value  of  a  criterion 
of  sensitivity  of  any  concentration  o*  with  respect  to  the 
same  constant  k  may  in  the  course  ofJthc  reaction  ohange  oon-' 
cidorably.  Hence  it  follows  that  for  a  comparison  of  tho 
values  of  criteria  of  sensitivity  for  various  o,  or  L,  it  is 
nccncso’v  to  stipulate  the  time  or  at  least  the^stage^of  the 
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reaction  at  which  the  comparison  is  being  made. 

As  could  have  been  expected,  for  one  moment  of  time  and 
under  the  same  conditions  of  the  reaction,  different  criteria 
of  sensitivity  have  different  numerical  values  but  of  the  same 
sign  (in  particular,  they  simultaneously  cross  the  x-axis)  (see 
Figures  i?a,  18-21). 

If  the  curves  corresponding  to  some  criterion  have  a  posi¬ 
tive  maximum  (or  negative  minimum),  then  the  curves  for  the 
other  criteria  may  have  the  same  characteristic  at  a  location 
shifted  some  vertical  and  horizontal  distance  away.  (See 
Figures  l?a3  10  and  11.)  The  criteria  of  each  of  the  concen¬ 
trations  of  CH20,  H20,  CH4,  Oo  and  CH~  with  respect  to  a  change 
in  various  constants  of  reaction  rateJk  at  various  stages  of 
the  reaction  are  shown  in  Figures  17a,  I7b,  18-20, 


Let  us  look  at  the  behavior,  for  example,  of  the  criteria 
of  sensitivity  of  the  above-mentioned  concentrations  at  a  time 
close  to  t~  =  2*10“3  sec.,  where  the  concentration  of  CH^O  reaches 
a  maximum  (see  Figure  17a) .  W*  note  that  the  concentration  of 
formaldehyde  at  this  stage  of  the  reaction  is  most  sensitive  to 
a  change  in  the  constants  of  reaction  rates  k/j,  and  ko  (which 
determine  the  rate  of  reaction  of  CH4  with  the  radicals  OH  and 
HO  )  and  in  the  constants  kr  and  kg.  The  sign  of  the  criteria 
of2sensifcivity  of  the  concentration  of  CHgO  to  a  change  in 
and  kg  is  positive,  while  the  sign  of  the  criteria  with  res¬ 
pect  to  a  change  in  and  is  negative.* 

* It  can  be  seen  in  Figure  17a  that  an  increase  in  fy,  (or  kg) 
of  about  10#  at  time  tffi  =  2*10“^  sec.  leads  to  an  increase0 
in  the  concentration  of  CH20  of  about  30/5 ;  at  the  same,  time 
an  increase  in  kQ  of  10,#  causes  a  decrease  in  the  concentra¬ 
tion  of  CM20  of  about  70#. 


At  the  same  time  the  sensitivity  of  the  concentration 
of  CH^O  to  a  change  in  the  constant  k,  (determining  the 
reaction  rate  of  CH20  with  02)  is  small  near  tm  =  2*10“3 
sec.  Tt  can  be  seen  in  Figure  17b  that  near  t_  =  2*10“3 
sec.  the  sensitivity  of  the  concentration  of  formaldehyde 
to  a  change  in  constants  k^,  k2,  kg,  k„,  k,.  (see  Table  7) 
is  relatively  small.  (Near  this  point 'of  time  several  of 
the  values  of  these  constants  change  sign.)  Near 
t  -  2* 10  -  sec.  the  concentration  of  Up 0  has  the  greatest 
sensitivity  to  a  change  in  kg  (see  Figure  18).  The  concen¬ 
tration  of  CHg  has  the  maximum  sensitivity,  with  a  positive 
sign  for  the  criterion,  to  a  change  in  constant  k and  a 
somewhat  smaller  absolute  magnitude  (but  with  a  negative 
sign  for  the  criterion)  of  sensitivity  to  a  change  in  kr- 
•'see  Figure  20).  It  is  possible  to  compare  the  sensitivities 
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not  reproducible 


of  the  concentrations  of  02  and  CIIi,  to  a  change  in  the  constants 

*<2,  311(1  ^8  ^see  Figure  19)  • 

It  is  also  possible  to  see  from  these  figures  the  extent 
to  which  the  relationship  of  the  values  of  the  criteria  of  sen¬ 
sitivity  of  one  concentration  with  respect  to  various  con¬ 
stants  changes  over  the  course  of  the  reaction.  In  this  instance 
there  are  several  cases : 

a)  The  qualitative  relationship  (less-more)  is  main¬ 
tained  over  the  i*hole  course  of  the  reaction,  but  the  quanti¬ 
tative  relationship  of  the  values  of  the  criteria  differs  in 
various  stages  of  the  reaction  (for  example,  it  can  be  seen  in 
Figure  17a  that  the  sensitivity  of  the  concentration  of  CH?0.to 
a  change  in  lq,  and  ko  is  greater  than  the  sensitivity  of  tnis 
concentration  to  a  change  in  kg,  for  the  whole  duration  of  the 
reaction). 

b)  The  qualitative  relationship  changes  in  the  course  of 
the  reaction.  For  example,  in  Figures  17a,  17b  it  can  be  seen 
that  this  is  the  case  for  the  sensitivity  of  CH20  with  respect 
to  a  change  in  the  constants  k^  and  k^;  k2  and  R^,  etc. 

Figures  17a,  19  and  20  show  for  purposes  of  comparison 
the  curves  illustrating  the  change  with  respect  to  time  of  cri¬ 
teria  (66)  and  (68);  while  Figure  21  shows  the  same  for  cri¬ 
teria  (68)  and  (69).  It  can  be  seen  from  these  figures  that,  as 
could  have  been  expected,  with  the  passage  of  time  these  curves 
behave  differently.  This  behavior  is  easily- explained.  For 
example,  the  difference  in  the  behavior  of  criteria  (68)  and 

in  Figure  21  is  explained  by  the  fact  that  criterion  (68) 
■determines  the  change  of  the  concentration  relative  to  its  cur¬ 
rent  value,  while  criterion  (69)  determines  the  change*  of  the 
concentration  relative  to  its  maximum  value. 

’.Jo  also  note  that  the  computations  that  have  been  car¬ 
ried.  out  enable  us  to  separate  out  those  constants  whose  changes 
have  the  least  effect  on  one  or  another  concentration  over  the 
whole  course  of  the  reaction.  Thus  it  can  be  seen  from  Figure 
17b  that  the  curves  of  dependence  on  time  of  the  criteria  for 
and  '<n ,  for  the  entire  reaction,  are  closer  than  other 
ctirv^r  tb  the  axis  of  time.  This  means  that  changes  in  con¬ 
stants  k~  and  k7  have  the  least  influence  on  the  concentration 

ryn  r**T  a  * 

Computations  showed  that  a  change  in  constant  k3  had 
little  influence  on  other  concentrations  as  well  (with  the 
exception  of  the  concentration  of  CH  00,  whose  magnitude  is, 
however,  small  itself).  This  circumstance  indicates  that  in 
the  reaction  process,  the  stage  that  is  determined  by  the  oon- 
:taru  k~  occurs  with  great  speed,  3uoh  a  oonolqsion  agrees 
with  tlw  results  obtained  by  Semenov  (souroe  T3ij  ),  From 
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the  point  of  view  of  the  mathematical  structure  of  equation  set 
(70),  the  small  influence  of  constant  k~  is  explained  by  the 
fact  that  this  constant  enters  the  equations  only  in  the  form  of 
the  derivative  of  ko  with  respect  to  c,,  and  the  latter  is  small 
over  the  whole  course  of  the  reaction. 

In  conclusion  we  shall  make  several  observations. 

1.  The  criteria  of  sensitivity  that  have  been  constructed 
must  permit  us  objectively  to  evaluate  the  role  of  one  or  another 
constants  of  reaction  rates  of  individual  stages  of  complex 
chemical  reactions.  At  the  same  time  we  find  which  reactions  are 
determining  (cruciaQ  for  the  process  in  its  various  stages. 
However,  we  do  not  yet  consider  the  interpretation  of  the  curves 
obtained  above  to  be  clear  enough. 

It  is  necessary  to  keep  in  mind  that  both  in  the  mechanism 
of  the  complex  reaction,  as  well  as  in  the  set  of  differential 
equations  corresponding  to  this  mechanism,  a  change  in  one  of 
the  constants  k  causes  in  the  system  as  a  whole  other  changes 
which  affect,  strictly  speaking,  all  of  the  concentrations  to 
one  extent  or  another  over  the  course  of  the  whole  reaction.  In 
determining  with  the  help  of  criteria  (65)-(6 9)  the  sensitivity 
of  various  concentrations  to  a  change  in  the  constants  k  .  we 
are  really  determining  the  sensitivity  of  the  whole  chemical 
process  to  that  stage  of  its  mechanism  that  is  characterized  by 
the  given  constant  k  .  At  the  same  time  it  can  not  be  for¬ 
gotten  that  because  of*  the  very  character  of  the  set  of  dif¬ 
ferential  equations  of  the  process,  in  any  such  criterion 
there  is  reflected  to  some  degree  or  other  the  changes  in  all 
of  the  concentrations;  that  is  to  say,  there  is  reflected  the 
indirect  influence  of  the  values  of  the  other  constants  kq. 

Prom  the  standpoint  of  cybernetics  the  reaction  mechariism4is 
a  system  whose  elements  are  mutually  interdetermined,  both  at 
any  instant  of  time  t,  and  also  in  the  system's  evolution  over 
t  ime , 


A  knowledge  of  the  magnitude  <ic./  “&k  is  essential  for 
the  proper  organization  of  the  experiment  for  determining  the 
constants  of  reaction  rates  k  and  for  evaluating  the  quality 
of  the  results  obtained. *  ^ 

•"■In  carrying  out  isothermal  reactions,  for  example,  the  terri- 
perature  in  determined  with  a  certain  amount  of  error  Axg. 

'Tier  the  consequent  error  in  the  determination  of  the  constants 
&k0  -  (dk„/dT)  A" .  Using  the  urreniun  equation  ’with  the 
energy  of  activation  3T,  we  find 
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TTencc  it  follows  that 
1. 


ACj/cj  may  be  quite  large  even  when 


The  criteria  of  sensitivity  may  be  used  not  only  in  the 
case  of  homogeneous  gaseous  reactions,  but  also  in  any  other 
types  of  chemical  processes  (heterophasic  processes,  processes 
occurring  in  condensed  phases,  solid  state  processes,  etc.). 
This  approach,  analogous  to  that  developed  above,  may  be  used 
for  the  study  of  the  influence  of  various  parameters  (for  exam¬ 
ple,  the  degree  to  which  the  surface  of  a  catalyst  is  occupied, 
the  length  of  a  molecular  chain,  etc.)  on  the  course  of  a 
chemical  reaction. 


2.  Tn  non-ieotherraal  reactions  it  is  possible  to  for¬ 
cibly  alter  the  temperature  of  this  or  that  stage  of  the  reac¬ 
tion  (it  can  be  done,  for  Instance,  in  the.  case  of  reactions 
occurring  in  a  c tat ionary  current).  Then  a  knowledge  of  the 
’’’ag'ii  tude  o**  Be.,/  3k_  will  help  in  the  evaluation  of  the  con- 

-  r* i * -  - -  - -  *-•-  ■  —  «  •  i  * 
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such  an  extraneous  thermal  action  on  the  coure-: 

'  complex  chonical  reaction.  A  similar  situation  arises 
rfVv*  it  is  necessary  to  take  into  account  the  influence  of 
1 1;-**  .;sllc  of  a  reaction  chamber  ,  additional  radiation 

outside  the  system  ,  or  additional  quantities  of  one  or 
another  reagents  introduced  into  the  area  of  reaction. 


*  - 


*:  •  .  'ft#*  "  •  • 

'he  approach  that  we  have  developed  may  be  extended  to-  the 
sf  non- isothermal  reactions.  For  this  .purpose  it  is  only 

IV 
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to  look  at  the  magnitudes  of 
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3c,/fckoa  and  3c,/ *>2n 
?  of  t>c,/3kg,  whe^e  kQq  his  the  pre-  q 
and  ic  the  energy  c reactivation, 
uf  equations  of  t)*pe  (71)  then  increases,  but  the 
the  equations  remain  the  same. 
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